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CN ; Abstract 

Consider a continuous time random walk in Z with independent 
and exponentially distributed jumps ±1. The model in this paper 
consists in an infinite number of such random walks starting from the 
complement of {— m, —m + 1, ...,m — l,m} at time —t, returning to 
the same starting positions at time t, and conditioned not to intersect. 
This yields a determinantal process, whose gap probabilities are given 
by the Fredholm determinant of a kernel. Thus this model consists 
of two groups of random walks, which are contained into two ellipses 
(■*") ' which, with the choice m ~ It to leading order, just touch: so we 

have a tacnode. We determine the new limit extended kernel under 
the scaling m = \2t + at 1 ' 3 J , where parameter a controls the strength 
of interaction between the two groups of random walkers. 
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1 Introduction 



In the past decade, systems of vicious random walks and non-intersecting 
Brownian motions have been investigated, and quantities such as the cor- 
relation functions [38J, the one-point distribution functions and limit pro- 
cesses under appropriate scaling limits have been studied. Non-intersecting 
Brownian motions arise in the study of random matrices J32J [33j |37] , and 
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space (and/or) time discrete versions in random tiling and growth mod- 
els I23H2HI2SH2Sl[3niSIlS23 - Most of these works use the mathematical frame- 
work shared by Brownian motions starting from a point, and either ending 
at the same point after a given time or the boundary condition is free (with 
possible extra boundary conditions like staying positive [36|l50]). 

Consider N non-intersecting Brownian bridges Xj(t) on R, leaving from 
at time r = —2N and forced to at time r = 2N. For large N, the 
mean density of Brownian paths has support, for each — 2N < t < 2N, on 
the interval (— y/AN 2 — r 2 , \/AN 2 — r 2 ). This means that on the macroscopic 
scale, where space and time unit are set equal to N, one sees a circle. Near 
its boundary, the density of Brownian paths is of order iV -1//3 , thus to see 
something non-trivial one needs to look in a space window of size N 1 ^ 3 and, 
by Brownian scalings, a time window of size N 2 / 3 . We call this the "Airy 
microscope" , since it holds 

lim P (all N- 1/3 (x t (2sN 2/3 ) - 2N) E E c - s 2 ) = P (A 2 (s) n E = 0) , 

(1.1) 

where A 2 is the so-called Airy 2 process. It has a universal character and 
was discovered in the context of the so-called multilayer PNG model [H] . 
The scaling (II. ip is equivalent to the customary iV _1 / 6 -GUE-edge rescaling 
along the circle for non-intersecting Brownian motions leaving from the origin 
at time t — and returning to the origin at time t — 1; this is done by an 
appropriate change of the variance of the Brownian motions. 

In the context of growth models, generalizations have been introduced 
with external sources [TTJ[25] • Its analogue in terms of Brownian motions is 
to require that a finite number of Brownian motions end up at some point 
aN. Then under the scaling in (II. ip . the limit process is a transition process 
from Airy 2 to Brownian Motion. For extensions to more general sources, 
see [9],[T8], while for the case that the top r Brownian motion end up at 2N, 
see [3] and [I]. 

A further known situation occurs when a fraction pN of the iV non- 
intersecting Brownian motions (leaving from the origin at time t — 0) end at 
time t — 1 at position aN and another fraction (1 — p)N at bN, with a < b. 
When iV — > oo, the mean density of Brownian particles has its support on 
one interval in the beginning and on two intervals near the end. Thus a 
bifurcation appears for some intermediate time To, where one interval splits 
into two intervals, creating a "heart-like" shape with a cusp at the origin. 
Near this cusp appears a new universal process, upon looking through the 
"Pearcey microscope" , where the space window is iV 1//4 and the time window 
is N 1 / 2 . The new process is is called the Pearcey process [17] and is 
independent of the values of a, b and p; see [§]. Once the bifurcation has 



taken place, the Brownian motions will eventually fluctuate like the Airy 2 
process near the edge, with a transition from the Pearcey to the Airy 2 process 
[2] . The Pearcey process has been obtained as the limit of discrete models, 

see mmang. 

The motivation of our work was to understand what happens when half 
of the non-intersecting Brownian motions start and end at a point, while the 
second half start and end at another point. When the two starting points 
are sufficiently far apart from each other, the mean density of particles will 
be confined to two separate circles, with Airy 2 processes appearing near the 
boundary, as described above. When the two starting points move away from 
each other at an appropriate rate proportional to N, the two circles will just 
touch, creating a tacnode. A new critical process appears by looking at 
the two sets of non-intersecting Brownian motions, which experience a brief 
meeting in the neighborhood of the tacnode, but looked at with the Airy 
scaling; we call it the tacnode process. Pictorially it can be thought of as 
two Airy 2 processes touching. 

In this paper we obtain an explicit formula for the kernel governing this 
tacnode process starting for non-intersecting continuous-time random walks, 
rather than non-intersecting Brownian motions. The same result is expected 
to hold for the Brownian motion case, since under the scaling the discrete 
nature of the random walks is lost and the random walks becomes Brownian 
paths. Our main result is the limiting kernel at the tacnode under appropriate 
scaling limit, stated in Theorem 12.21 Before taking the limit, the kernel is 
given by Theorem 12.11 The model is to let two groups of non-intersecting 
random walks with jumps ±1, rate 1 and 2m + 1 integers apart evolve during 
a total time or order m, with space-time rescaled a la Airy, namely x ~ £m 1//3 
and t ~ sm 2//3 as suggested by formula (II. ip . 

There is an important difference with respect to the previous two cases: 
here we have a one-parameter family of processes, which is obtained by modu- 
lating the end points distance between the two sets of Brownian motions over 
distance of order N 1 ^ 3 . For the Pearcey processes (and the Airy 2 process), 
geometric changes of this type only have the effect to modify the position 
(and orientation) of the cusp, but the underlying Pearcey process remains 
unchanged. In the literature there is another known situation with a process 
in a tacnode-like geometry [H] . which however differs from the present one. 

This problem can be approached using multiple orthogonal polynomi- 
als [20] and a Riemann-Hilbert to this problem is analyzed in a paper [21] 
(which meanwhile appeared in the arXiv). In a forthcoming paper [30] Jo- 
hansson uses a different approach leading to a different kernel for the Brow- 
nian motion problem, which is expected to be the same as the kernel in 
our paper. Adler, Johansson and van Moerbeke have then considered two 



partially overlapping Aztec diamonds and found the same kernel [5]. 

Outline 

In Section [2] we define the model and state the two main results. In Section [31 
Theorem 13.11 we derive the finite time result for r = 0, which is reshaped 
in Section [4] as a preparation to carrying out the large time limit. Before 
actually doing this, we indicate in Section [5] how to introduce the time, 
leading to the finite multi-time kernel in Theorem 15.41 an extension of the 
kernel appearing in Proposition 14.11 In Section [HI we take the limit of the 
multi-time kernel, leading to the proof of the first formula of Theorem 12.21 
In Section we sketch the proof of the double integral representation of the 
kernel, the second formula of Theorem 12.21 using the steep descent analysis. 
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2 Model and results 

Consider a continuous time random walk in Z with jumps ±1, occuring 
independently with rate 1, i.e., the waiting times of the up- and down-jumps 
are independent and exponentially distributed with mean 1. The transition 
probability Pt(x, y) of going from x to y during a time interval of length t is 
given by 

p t (x,3/)=e-%_ 1/ |(2t), (2.1) 

where I n is the modified Bessel function of degree n (see pQ). 

Consider now an infinite number of continuous time random walks start- 
ing from {. . . , — m — 2, — m — l}U{m+l, m+2, . . .} at time r = — t, returning 
to the starting positions at time r = t, and conditioned not to intersect, see 
Figured! Denote Xk(r) the position of the walk that starts and ends at posi- 
tion k. Then, the point process ijonZ (described by the little white circles 
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Figure 1: The lines are the non- intersecting walks x. The white dots are the 
support of the point process fj. 



in Figure [T]) defined by 



fj(x) 



E 

fc€Z\{-m,. 



Vx,x k (0)i 



(2.2) 



,m\ 



with 5 the Kronecker-delta, is determinants!, i.e., there exists a kernel K m 
such that the /c-point correlation function p^ is given by p( k \yi, . . . ,yk) = 
det(K m (yi,yj))i<i t j<k- One of the interesting quantities is the gap probability 
of a set E, which is given by P(jj(1e) = 0), i.e., the probability that none of 
the random walks are in E at time r = 0. For a determinantal point process 
the gap probability is given by the Fredholm determinant of the associated 
kernel K m projected onto E. For more informations on determinantal point 
processes, see [EMIEHIS! ■ 

The determinantal structure still holds if we consider the point process 
on a set of time-slices instead of a single time r = 0. This means that 
given times T\ < r 2 < . . . < r p in the interval (—t,t), the point process on 
{n, . . . , t p } x Z defined by 



fj(r,x) 



E E 

r=l fceZ\{-m,. 



8(t,x 



(T,x),{T r ,X k {T r ))l 



(2.3) 



.,m} 



is determinantal. That is, the space-time correlation functions are given by 
the determinant of an extended kernel, which we denote by K^(ti, x\\ t 2 , x 2 ), 
where tj G {ti, . . . , r p } and x, G Z. 

It is more convenient to first study the dual or complementary process 
x(t). The dual proceeds along the gaps of5t(r). In this instance, the dual x(r) 
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Figure 2: The solid lines are the non-intersecting walks x, the dual process 
of x of Figure HJ The black dots are the support of the point process r\. 




Figure 3: Superposition of Figure [Hand Figure [2] 



of x(r) is described by n = 2m + 1 (m e N) non-intersecting continuous-time 
random walks, starting from — m, —m + 1, . . . , m at time r = —t, returning 
to the starting positions at time r = t; see Figure El and Figure [3] for the 
superposition of the trajectories of x(r) and x(r). 

In particular, the dual process x(r) at r = is given by the the little black 
circles in Figure [2J The probability measure at time r = is obtained by 
the Karlin-McGregor formula [21], and thus it is a determinantal process for 
a kernel K m . Finally, the complementation principle by Borodin, Olshanski, 
and Okounkov (see Appendix of [T7j) tells us that, if the kernel K m governs 
the process x(r), then the kernel K m = 1 — K m describes the dual process 
Mr). 



Theorem 2.1. The determinantal point process ^(r, x) on {r l5 . . . , r p } x R, 
Ti G (— t, t), defined by the two groups of non-intersecting walkers, starting 
and ending 2m + 1 apart, at times —t and t respectively, has gap probabilities 
on any compact set E C {ti, . . . , r p } x R gro en by 

P(fj(l E ) = 0) = det(l - K^ xt ) L2(i ,), (2.4) 

where the kernel K^f zs gzi> en by 

e 2t2 ~ 

—K e * t {t u x 1 ;t 2 ,x 2 ) = -l [t2<tl] I\ Xl „ X2 \(2(t 2 - t x )) 



(27ri) 2 Jp 7 roz e t(w-w-^) e -t 2 { w + w -^) zXl - m z-w 

Vm r dw I dz 6 ^'"'^ e ~ tl{z+Z ~ 1] wX2+m H 2m+1 (z)H 2m+1 (w- 1 ) 



(27ri) 2 Jp J r e t(z-z i) e -t 2 (»«+«'- :L ) ^i+m+1 

i^i^lT^T f d z^--——T—H 2m+ i{z- 1 )H 2m+l (z), 



w — z 



27ri./ r " z xl ~ X2+1 



i) 



(2.5) 

TOi/i Kn := l/(-^2m+i(0)i^2m+2(0)). T/ie function H n is itself the Fredholm 
determinant on £ 2 ({n, n + 1, . . .}) 

H n (z~ r ) := det(l - ^(z- 1 ))^^,^!,...}) (2.6) 

of the kernel 

(—l)k+£ r r u e 1 u — z e 2t{ - u ~ u ~ 1 ^ 
Kiz- 1 )^ := ^— %- f dui dv^— 27(^TT> ( 2 - 7 ) 

where To is any anticlockwise simple loop enclosing and similarly Tq }U en- 
circles the poles at and u (but not z)j. 

The extended kernel, governing the process fj(r, x), is given in terms of the 
the kernel K m (xi,x 2 ) = K^ t (0, X\\ O,:^), governing the distribution fj(0,x), 
by 

K% t {t 1 ,x 1 ;t 2 ,x 2 ) = -t[t 2 < tl ] (e^ 2 - tl)n ) (x u x 2 ) + (e- hn K m e t2H ){xx,x 2 ), 

(2.8) 
where Ti is the discrete Laplacian 

(Hf)(x) = f(x + 1) + f(x - 1) - 2/(x). (2.9) 



1 For any set of points S, the notation JL dzf(z) means that the integration path goes 
anticlockwise around the points in S but does not include any other poles of /. 



Remark that the transition probability of (12. ip . defined for t > 0, can be 
written as p t (x,y) = e m t(x,y) =: e m (x,y). Here, 1 denotes the identity 
operator on Z, i.e., l(x, y) = 1 if x = y and t(x, y) — if x ^ y. 

The formula for the kernel K. m (ti,xi, £ 2 , X2) (that is at ti = 0) in The- 
orem 12.11 will be established in Section |3j whereas the one for K^* will be 
shown in Section [51 In Sections H] and [51 it will be shown that both kernels 
K m (x,y) and K^ t (ti,x 2 ; t 2 ,x 2 ) have a representation, whose constituents 
can be expressed in terms of Bessel functions; see the expression (14.141) and 
the time-dependent kernel (I5.26p . derived from (14.141) . via the recipe ()2.8)) . 
Also, note that the kernel K(z~ v ) is a rank-one perturbation of the kernel 
K(0), whose Fredholm determinant 

H n (0) = det(l - K(0)) £2({nin+1 ,... }) (2.10) 

is the distribution of the longest increasing subsequence of a random permu- 
tation in the Poissonized version, or, equivalently, it is the height function in 
the polynuclear growth (PNG) model [T014T] . In the scaling limit, considered 
in Section El H n (Q) will converge to the Tracy- Widom distribution F 2 . 

To study the limiting behavior, when m, t — >• 00, consider first the system 
of non-intersecting random walks starting at time —t and ending at positions 
{. . . , — m — 2, — m — 1} at time t. This is, up to a shift by m+ 1, the multi- 
layer PNG model studied by Prahofer and Spohn in [41 J. Their work shows 
that the top random walk at time r = has fluctuations around x = —m + 2t 
of order t 1 ' 3 . By symmetry, if one considers only the non- intersecting random 
walks starting and ending at position {m + 1, m + 2, . . .}, the bottom random 
walk at time r = fluctuates around x = m — 2t also in the spatial scale 

The top and bottom random walks interact if the proportion of deleted 
configurations, due to interaction, is non-zero. This happens when m — 2t to 
leading order in t. The first scaling where interaction is relevant is given by 
m = 2t + at 1 / 3 . The parameter a modulates the strength of interaction of the 
two sets of non- intersecting random walks. In the extreme cases a — > 00, we 
clearly (by a simple probabilistic argument) go back to the situation of two 
independents PNG models, thus the top of the lower walks and the bottom 
of the upper walks are governed by the Airy 2 process [41 J. On the other 
hand, when a — > —00, one expects to see a point process governed by the 
sine kernel or the Pearcey process. Moreover, locally the paths will looks like 
random walks, so the exponents in the scaling for time and space are in a 
ratio 2:1. Thus, we set the scaling^ 

m = 2t + at 1/3 , x t = & 1/3 , U = Sit 2/? % i = 1, 2. (2.11) 



2 We do not write explicitly the integer parts, since in the t — > 00 it is irrelevant. 



Also note that for each time —t<r<t, the density of particles has its 
support on two semi-infinite intervals, whose boundary, as a function of r, 
describes two curves, which at r = form a tacnode. The purpose of The- 
orem 12.21 is to describe the fluctuations of the random walks in the t — > oo 
limit in the neighborhood of (x, r) = (0, 0), but in the new space-time scale, 
given by (12. lip . 

In order to state the second main result, define the standard Airy kernel, 



#Ai(6,6) == / dXAifa + A)Ai(6 + A). (2.12) 

Jo 

and the function Q(k), already appearing in [48J, 

Q{k) := [(1 - XaKMXar'xaAi]^), with a := 2 2 / 3 a, (2.13) 

and where Xa( x ) = l{x>a]- We further set 

Ai (s) (0 := e Cs+ f s3 Ai(e + s 2 ), (2.14) 

which equals to the standard Airy function Ai(£), when s — 0, and define 
the functions 

/•oo /»oo 

A(s,0 :=Ai (s) (a-0+ / dK da Q(k) Ai(k + a) Ai {s) (2 1/3 a + a -£), 

Ja JO 

poo 

B(s, := / dK Q( K )Ai^(2 1 / 3 /€ - a + £), 
J a 

(2.15) 

and 

/■oo r 

C(s,0:=2- 1/3 / d«Q(K) Ai( 2 " /3 ^( K + 2- 1 / 3 

/*00 /'OO 

+ / d\Q(X) daAi(a + \)Ai {2 ~ 2/3s \a + K + 2- 1/3 + (f ** -£), 

(2.16) 

where with (£ <H- — £) we mean the same expression with £ replaced by — £. 
Finally, we define two Laplace transforms V(u) and Q(u): 



Q{u) :- 



dn Q(n)e 



Ku2 x l' i 



V(u) := — / due 
'o 



-Ku2 x / 3 



(2.17) 



dfi Q(/i)Ai(/i + k). 



Theorem 2.2. Near the tacnode appears a new determinantal process on 
{si, . . . , s p } x R, the tacnode process T , whose gap probabilities on any com- 
pact set E C {si, . . . , s p } x R are groen fry 

P(T(1 B ) = 0) = det(l - /C cxt ) L2(£) . (2.18) 

The kernel /C ext zs ine limit of K^f under the scaling 112. 11\) . 

( — \) x 2p4t2 _ 

/C cxt (si,£i;*2,6) := lim t^K^s^Sa), (2.19) 

t— >oo ( — M xi e * 

where the convergence is uniform for £i,£ 2 one? si,s 2 in bounded sets. The 
kernel /C ext nas the following representations: 

/C ext (si,£i;* 2 ,6) = - /A 1[S ; <S1] ex P f- 5 1 " 6 ^ ) + C(*i -*.,&- 6) 

+ f d-f(A{ Sl ,^ - 7)^(-s 2 , 6 - 7) + ^(si, -6 - 7)^(-s 2 , -& - 7) 
-^1,6 - 7M-S2, 6 - 7) - A{ Sl , -6 - 7)^(-«2, -6 - 7) 
-B(8i, & - l)A{-s 2 , 6 - 7) - B(si, -6 - -y)A{-82, -£2 - 7) 

- /" d 7 r j B( 5l ,ei- 7 )S(-S2,e2-7)+B(si,-ei-7)B(-«2,-6-7)), 

(2.20) 
as we// as 

£ si, ft; s 2 , 6 = 7= . exp -— + C(si - s 2 , 6 - 6) 

J Ait si - s 2 V 4(si-s 2 )/ 



• 3 2 



1 /",/", e-'^e 81 " /e ?1 " e^ iu \ (1 - P(n))(l - P(-n)) 



(27ri) 2 7 5+iK i_ 5+iM e 4- ff , e^ 2 \e?* v e-&>J u-v 

1 A r elf"™ e siu2 /e« lU e- ?lU \ (1 - P(n))Q(-t>) 

(2vri) 2 7 25+ « X + « e -£-<™ e^ 2 (^ + 1^) u~^v 



" 3 „2 



CL Hi i Uj U o IT" I T I 



(27ri) 2 J-5+iR J-25+m e 4-« e s ^ 2 Ve«^ e-&V u-v 



■• 3 ,.2 



1 [ r [ 7 e -^-^ e SlU ( e^ u e~^ u \ Q(u)Q(-v) 

du dv 5 5- — : h 



(2vri) 2 J_ s+m J 5+m e -$-ov e^ 2 \e^ e~& J u-v 

(2.21) 

The form ( I2.20p of the limiting extended kernel in Theorem 12.21 will be 
shown in Section [61 whereas a sketch of the proof of its double integral 
representation (I2.2ip will be given in Section [7J 

10 



In the preprint [21] the analogue problem for Brownian Motion will be 
analyzed with the Riemann-Hilbert approach to multiple orthogonal polyno- 
mials. It would be interesting to see how to relate the two formulas (which 
we expect to be the equivalent). 

3 Finite system at r = 

In this section we will prove Theorem 12. II for t\ — t^ — 0. Consider a contin- 
uous time random walk in Z with jumps ±1, which occur independently with 
rate 1, i.e., the waiting times of the up- and down-jumps are independent and 
exponentially distributed with mean 1. Thus, the number of up-jumps (and 
similarly down-jumps) during the time interval [0,t] is Poisson distributed, 

_ t k 
P(k up-jumps during [0, £]) = e — . (3.1) 

As will be shown, the transition probability pt(x,y) of going from x to y 
during a time interval of length t is given by 

p t (a;,3/) = e-%_ v |(2t), (3.2) 

where I n is the modified Bessel function of degree n (see PQ). To prove (13. 2p . 
first notice that by symmetry, it is enough to consider y — x > 0. To go 
from x to y, the process must perform k steps to the left and k + y — x steps 
to the right. Since the moment, at which the left or right steps occur, is 
independent of whether it is a left or a right step, on may assume the process 
doing first k steps to the left and then k + y — x steps to the right. By the 
strong Markov property of the random walk and the independence of the 
jumps, 

Pt (x, y) = }_,P^ k down . steps j during time tj 

The modified Bessel function has the following expressions (for n G Z) 

I n (2t) = — J> - e ^V^ y , m , (3.4) 



with S 1 = {zeC\\z\ = 1}. 

11 



l<i,j<2m+l L^Vi")"" ' - J I Jl<ij<2m+1 

2 



Consider now n = 2m+l (m e N) continuous time random walks starting 
from — m, — m + 1, . . . , m at time r = — £, returning at the starting positions 
at time r — t, and conditioned not to intersect. Denote by Xkij) the position 
at time r of the random walk which started from m + 1 — k (i.e., the fcth 
highest one), see Figure [2] for an illustration with m — 2. 

The probability at time r = is easily obtained by the Karlin-McGregor 
formula J3TJ, namely 

(2m+l 2m+l 

f| {x*(0) = y fc }| f| ^ Xfc W = Xfc ( _t ) = m + 1 - A;} 
fe=i fc=i 

= const x det [p t (m + 1 - i, %•)] 1 < iJ < 2m+1 det [pt(jfa, m + l-j 

= const x (det [J K+j _ 1 _ m (2t)] 1 <.^ 2m+1 

(3.5) 
It is well known by [13] that the process above 

x(r) := {x fc (r),l < A; < 2m + 1}, re[-t,t], (3.6) 

with a measure of this form gives rise to a determinantal point process (ran- 
dom point measure) 

2m+l 

fc=i 

with a certain kernel K m (x,?/), to be computed in Theorem 13. II 

Instead of the process x(r), we shall analyze its complementary (dual) 
process, which we denote by 

x(r) = {x k (r),ke Z\{1 < fc < 2m + 1}}, tG[-t,t]. (3.8) 

If x denotes the trajectories of the 2m + 1 particles, then let x denote the 
trajectories of the holes, obtained by the particle-hole transformation, see 
Figures [T] and El 

The reason for starting with the process x is that the Karlin-McGregor 
formula applies to a finite number of paths, while x has an infinite number of 
paths. By the complementation principle in the Appendix of [17], the dual 
point process at r = 0, 

f l = J2 5 x k (o), (3.9) 

k 

is also determinantal with correlation kernel 

K m (x, y) = 5 X)V - K m {x, y). (3.10) 

First of all, we compute the kernel K m (x, y) in a form which will be suitable 
for asymptotic analysis. 

12 



Theorem 3.1. The point processes r\ and fj, defined in (ff.7|) and Ii3.9\) . are 



determinantal with correlation kernel K m and K m given below. Thus, for 
any finite subset Eel*, the gap probability of E is given by 

P { V {1 E ) = 0) = det(l - K m ) m , P (fj(t E ) = 0) = det(l - K m ) P{E) , 

(3 ' 11} 
with kernels K m (x, y) and K m (x,w) ; invariants under the involution 

(x,y) -H- (— y, — x), namely 

^m(x,y) = — £& fcjv, ^—^ — __ 

+ K> f d f d /^ ***» H 2m+1 ( Z )H 2m+l{w -i) 



(2vri) 2 / ro J Tqw e t(z-z ^) z x+m+i w-z 

+ 77^9 dz-——H 2m+1 (z- 1 )H 2m+1 (z), 
2m J To z x y+i 



(3.12) 



and 



at , , v m r r e'^" 1 ) w y- m ~ l H 2m+1 {w)H 2m+1 { z - 1 ) 

K m (x,y) =- (p dz cb dw- 



(27ri) 2 J ro J Toz e t{w-w *) z x-m z _ w 

Vm f d W I dz^^^ ^^ H 2rn + l(z)H 2m+1 (w~ l ) 



(2vri) 2 J Fo J Tow e t(z-z ^) z x+m+l w-z 

-M&vlT^f dz x - v+1 H 2m+i{z~ 1 )H 2m+ i{z), 

Z7ri j ro Z y ^ 

(3.13) 
where V m = l/(H 2m+ i(0)H 2m+2 (0)). The function H n itself is a Fredholm 
determinant on £ 2 ({n, n + 1, . . .}) 

H n (z~ l ) := det(l - K{z- l )) P{{n ^... }) (3.14) 

of the kernel 

(_l)fe+^ f f u e 1 u _2e 2i ( M_1 > 

* ( *~ 1)w := 727IF 1 rfM l, u dv ^^—z#*^*> ' (3 ' 15) 

where T is any anticlockwise simple loop enclosing and similarly F 0jU en- 
circles and m on/y (hence not z). 



3 As it should from the geometry of the problem! The involution interchanges the two 
double integrals in (|3.12p . as is seen from renaming w ■<-> z in the second double integral; 
also the third term, the single integral, only depends on \x — y\, as is seen from z — > z~ x . 
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Proof of Theorem \3.1[ 

Step 1: Computing the kernel K m (x,y), from the Karlin-McGregor for- 
mula i3.5\h It is well known by [13] that a measure of the form (13.51) implies 
that the point process (random point measure) r], as in (J3.7)) . is determinantal 
with correlation kernel 

2m+l 

K m (x,y) = y^ (p k (y)[A~ l ] kie <p e (x), x,yeZ, (3.16) 

k,l=l 

where 

<Pk{x) = 4+*-i-m( 2 *) (3.17) 

and A is the (2m + 1) x (2m + 1) matrix with entries 

[A] k ,t = (<Pk, <fe) = ^2 ¥k(x)p e (x). (3.18) 

Using (13.171) . the entries of the (2m + 1) x (2m + 1) matrix A, as in (I3.18p . 
are given by 



x'GZ x>0 x<0 



A k ,t = ^2ipk{x)<pt{x) = y]fk(x)(pt(x) +y^ j ip k (x)(p e ( 

xGZ 

^ O.TT 



X>0 



pt(z+Z ^-tfw+W 1 ) 1 

dz (b dw- 



Wfho L — ^ — T^r^ ( 3 - 19 ) 



Jl{z+z 1 )„t(w+w 1 ) i 



+ ^ (2vri) 2 % Q dZ f To dW z k w* (zwy-™ ' 

In the first integrals, we deform the paths to \z\ = 1 and \w\ — R > 1. Then 
we take the sum inside the integrals and use J2 x>0 (zw)~ x = wz/(wz — 1). 
Similarly, in the second integrals, we deform the paths as \z\ = 1 and 
\w\ = 1/R < 1 and use Yu X <o( zw )~ X = ~ wz /( wz ~ !)• This leads to 



Au o = — ■ ® dz (b dw 



e t(z+z 1 ) e t(w+w 1 ) wz 



(27Ti) 2 /| z |=i J\w\=R z k-m w £-m wz _ I 

If f e t(z+z- x ) e t(w+w- v ) wz 

— <£ dz (b dw ; -. (3.20) 

(27ri) 2 j^| =1 J H=1/R z k- mw e-m wz _ 1 

dz .,,. = h-t{4t), 



4111 J\z\ = l Z 

since for any value of z, the two integrals differ only by the residuqj at 
w = 1/z. However, doing the asymptotics of the kernel K. m (x,y) with this 
choice of basis and thus with this A -1 seems to be hopeless. 

4 This residue argument will reappear later in (I3.40J) . 
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Step 2: Changing the basis cp k (->■ tf; k , such that A *— y 1 in the kernel 
~K m (x 1 y), i.e., so that K m (x,y) = Ylk™i ^*( a: )'0fcG/)- Replace the basis 
{<Pk{x))k=i,...,2m+i with an orthonormal basis ('0fe(^))fe=i,...,2m+i with respect 
to the measure dp t (z) := ^^e^ z+z ' (generating the same vector space, i.e., 
det((^fe(^))i<fcj<n = const x det(ipk(xj))i<k,j<n so that the measure (13. 5p 
has the same form but with A = 1). More precisely, we shall search for 
polynomials P k of degree k such that 

dp t (z 



ipk{x) 



i z* 



-Pt-iiz- 1 ] 



dpt(w)w 



■^Pk-^w), l<k<2m+l, 



IS 1 * JS 1 

satisfies, using the same argument as in (13.20 p . 

h,i = (il>kM = ^2f d Pt( z ) f dp t {w)(zw) x ~ m P k - 1 (z)P t -i{w) 

dp 2t {z)P k -i{z)Pe-i{z- 1 ) =: ((P k ^,P e ^)), 



(3.21] 



(3.22) 



,S'i 



thus defining a new inner-product (( , )) on the circle S 1 = {z £ C \ \z\ = 1} . 
So it suffices to find an orthonormal basis of polynomials on the circle for 
the weight dpitiz). A classical expression for the polynomial Pk(z) is (see 

e.g. my- 



Pkiz) 



•^/det m k ■ det m fc+1 



det 



/ 1 \ 

\Pi,j] 0<i<k Z 

o<i<fc-i 



(3.23) 



where m k = [pi,j}o<i,j<k-i and 



Z\7P) 



\ z k } 

dp 2t (z)z i - j = I i Jti) 



(3.24) 



Hence the P k (z) are polynomials of z with real coefficients. Orthonormal 
polynomials on the circle satisfy a Christoffel-Darboux-type formula, due to 
Szego; see [4"3] . Namely, with the notation P*(z) = z n P(z~ 1 ) and further 
using the reality of the coefficients, one obtains for z,w <E S 1 , 

n— 1 ra—1 



n(*)K(v) 



PJz)PJw) 



i=o 



1 — zw 



z n P n (z- l )w n P n (w- 1 ) - P n (z)P n (w) 
1 — w/z 

Z- n P n {z)w n P ri { W - 1 ) - Pnjz-^Pnjw) 

1 — w/z 
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(3.25) 



Step 3: Expressing the polynomials P n (z) in terms of the Fredholm deter- 
minant H^z' 1 ), as in ^S.lJjty - In order to do this, one first introduces the 
inner-product 

(f:9ks :=^i -fiuMu-^T^-^, (3.26) 

11X1 J s i u 

upon setting t := (ti, t 2 , . . .) G C°° and s := (s\, S2, ■ ■ .) G C°°. It was shown 
in J7J[8] (see also the lecture notes [51]) that the functions^ 

rf)(t, S ,):^^|±M= 

A/r n (t,s)r n+ i(t,s) 



(2)n , . .,- r n (t,s + [z ']) 



pW(t,s; *):=*" 



- if (3.27) 

r n (t,s)x„ + i(t,s) 



are bi-orthonormal polynomials with regard to the inner-product (I3.26p . In 
the formulae above, the r n (t,s) are 2-Toda r-functions and are defined as a 
Toeplitz determinant, which is also expressible as a Fredholm determinant 
of the kernel ( I3.29P below, using the Borodin-Okounkov identity [T6] . We 
obtain 



r n (t,s) := det 



27ri ./ s i m 



^V-^ e Er =1 fe« j -^- j ) 



KJW<n 



= Z(t,s)det(l-K(t,s)) £2({njn+W) , Z(t,s):=e-^ii*»-J, 

(3.28) 
where the kernel K(t, s) is given by 

Iff u l 1 e E ^ (i ^ J+ ^ v) 

K(t, b) m := ^ f rQ du ^ efo— —-^^-^-— . (3.29) 

The coefficients t,-, Sj have to be such that the expression YlT=i^j u ^ ~ SjU~i) 
appearing in the exponent of (13. 28p is analytic in the annulus p < \z\ < p^ 1 
for < p < 1. Then, the Borodin-Okounkov identity (I3.28J) gives a kernel 
K(t,s), with contours given by \u\ = \v | _1 = p', with < p < p' < 1. Using 
Cauchy's Theorem the contours may be deformed to any circle of radius 

< p < 1. Then, using Y1T=\{ V I Z Y 1 3 = ~~ m (l ~~ v / z ) (f° r \ v / z \ < 1) we 
obtain 

If f u l 1 1 - - e E ^i (tjt '~ J+s ^ J) 

K(t, s + [,-])„ = p-jj f^du £ ^——^-^^—^ 

(3.30) 



,2 „,3 



5 For a G C, one defines [a| = ( a, ^-, ^-, . . . ) G 
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and 

Z(t,s+ [z- 1 ]) = e-E^iit^-^-Vi) = Z(t,s)e-^i^ J '. (3.31) 

We now specialize all this to the locus 

On this locus, one checks that Z(t, s)| = e 4 ' , that K(t, s) and its translate, 
restricted to the locus C, are closely related to the kernel K(z~ l ) defined 

K(t,s)L c = j ir(0), 

V Jlc K ' (3.33) 

K(t,s+[^ 1 ])| £ ^ nj ^(^ 1 ), 

and that the restriction of r n (t,s) to £ leads to the Fredholm determinant 
H n {z~ l ) as defined in (13TT4J) : 



r„(t,s)| £ = fl n (0)Z(t,s)| jC = e 4t # n (0), 
r n (t,s+ [^ 1 ])| £ = ^"V^CM^ = Hniz-^e 4 *-*'*. 

Moreover, the inner-product (f,g)t,s defined in (13.261) reduces to the inner- 
product ((f,g)) defined in ([3T?2D . 

(f,9)t,s\c = ^i -e^+^fiuMu- 1 ) = ((f,g)). (3.35) 

"~ Z7T1 J s i U 

It follows that the bi-orthogonal functions for (f(z),g(z)) t)S restricted to the 
locus L coincide with the orthonormal polynomials defined by ( I3.22p . which 
by ([227D, (EE3D and (13331 yields: 

r n p -2t/zjT (~-\\ 

P n {z)=ti\t^z)\ c =ti\t,s-z)\ c = Ze Hn[Z , (3.36) 

^# n (0)# n+ i(0) 

where 

Fn^" 1 ) = det(l - K( 2; - 1 )),2 ({nin+li ... }) , (3.37) 

with the kernel K(z _1 ) as in ( 13.1 5ft ; this follows from ( I3.33p . The fact that 
the ph. and p n are equal on the locus £ is a consequence of the symmetry 



6 With A = B we mean that the two kernels A and B are conjugate kernels. In the 
present case, the conjugation factor is (— l) k ~ e . We remind that two conjugate kernels 
define the same determinantal point process. 
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of the inner-product (( , )), as in (13.221) . However, one easily verifies it with 
the above formulae. The equivalence of the Fredholm determinant parts is 
evident only after the change of variable v — > 1/u and u — > 1/v. Then, the 
kernel obtained for p n is the transpose of the one for p n . 

Step 4 : Expressing the kernel K m (x,y) as l[3.1ty) . Using this new basis ip^ 
as in f)3.2ip . and using the Christoffel-Darboux formula (13.25 p . the kernel 
K m (x,y) becomes by Step 2 (recall that n = 2m + 1): 



n r r „..y—m ™ r 

K m (x,y) = y2^k{x)^ k {y) = <p dpt(z) f dptiwj—^y^Pkiz^jPk 



' ' '-\ ( ir 
fc=l " J ' Jr °- z k=0 

= I dp t (z) I dp t (w) f "^ 1 ((-)" P n (z)P n (w- 1 ) - P n (z- l )P n (w) 
Jt Jy 0z z z — w \\z / 

(3.38) 

Note that the w-integrand in the double integral = has no pole at w — z, 

enabling one to deform the w-contour so as to include z G S 1 ; this has the 

advantage that the double integral of the difference can be written as the 

difference of two double integrals, each of them being finite. 

Inserting (I3.36P into (I3.38J) and setting V m = l/(H 2m+1 (0)H 2m+2 (0)) we 

get 

^ f ^ Vm I , / , e^-^ 1 ) w V-m-l H ( W )ff 2m+1 (^) 

K m (x, y) = ——r <b dz <p dw it— — r; — 

{2m) z J Vo J Toz e t{ - w w > z x m z — w 

V m f , f , e« w - w _1 ) W y +m H 2m+1 (z)H 2m+1 (w- 1 ) 



-, ^ dz w dw , . . 
(27ri) 2 J ro J Foz e t(z-z-i) z x+m+i z _ w 

(3.39) 
The expression in ( I3.12p is finally obtained by noticing that 

1 I j I j F(z,w) 1 1,1, F(z,w) J dz ^ , 

(M 5 t dz t dw ^T = J^f 1 *" t , dz ^T + 1 ^> *>• 

(3.40) 
proving formula (I3.12p . 

Step 5: Expressing the kernel K m (x,y) as A3.13\) . First of all, by (13. 36 p . we 
have 

H^z- 1 ) = P n (z)e 2t / z z- n ^H n (0)H n+1 (0). (3.41) 

Thus (with n = 2m + 1), the last term of ( I3.12p is given by 
l '" f <l: -H 2m+1 (z-')H 2m+1 (z) = -L I -^—e^+^P n (z)P n (z- x ). 



27ri J Fo zx-v+i""'^~ >-™r^~> 2m J sl Z-V+ 1 

(3.42) 
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In particular, 


at x = 


- y we have 










(|3.42|) 

» x— y 


= ((Pn, 


Pn 


and thus 












dz 

z x-y+l 


H 2m +l{z )H 2rn +l(z) 





1, (3.43) 



v m r i , (3 - 44) 

= 5 x>y + (l-5 x>y )— f dz———H 2m+l (z )H 2m+l (z). 

So, ¥L m (x,y) = 5 x>y — ~K m (x,y) = K^f(0, x\\ 0, x 2 ) of ( 12. 5p . thus establishing 
Theorem 13.11 This also ends the proof of Theorem 12.11 for t% = t 2 = 0. □ 



4 Reshaping, motivation and Bessel repre- 
sentation 

In this section we first reshape the kernel (12.51) of Theorem 12. II for t\ = t 2 = 0, 
to make it adequate for asymptotic analysis. Secondly, we rewrite all the 
terms using Bessel functions and the Bessel kernel. This will allow us to use 
known asymptotics for Bessel functions and kernel, without the need for new 
asymptotic analysis. 

4.1 Reshaping 

Note that the kernel K{z'~ 1 ) 1 defined in (J3.15J1 . with \u\ < \v\ < |z|, namely 
(_i\k+t r r 7 / i qi _ ? P 2t(u- U - 1 ) 

^ -1 )M := 7^V t du t dvJ m \t(v-^y ( 41 ) 

{2iti) z J To J Tou v k+i v — uv — z e zt( - v v > 

is a rank-one perturbation 

K{z-%t = K(0) k/ + h k (z~ l ) g e , (4.2) 

of the symmetridj kernel 

K(0) k , e = ij-L- I dui dv— —^-, (4.3) 

(2vri) 2 J Fo J Tou v k+l v-ue M{ - v v > 



upon using the identity 

1 u — z 



v — uv — z v — u v — z 



(4.4) 



7 as is seen by replacing u <—> 1/u, v n- 1/v. 
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where (remember \v\ < \z\ in the first integration below) 

-1 f dv e-^-v-^) 



hJz- 1 ) = — 



2vri /r„ (-v) k+1 







v — z 



dv e~ 2t ^ v ~ v ~^ e-m*-*- 1 ) 



2vri 7 ro 2 (-v) k+1 v-z (-z) k+1 



(4.5) 



--:h k {z-') + 



e 



-2t(z-z- 1 ) 



[-z) k+1 
and 



91 = 7^1 dui-uYe^-^. 
2vri J ro 



(4.6) 



In (I4.5p . one has replaced the integration about a small circle around by 
an integration about a contour containing z as well; this is done in order to 
be able to expand, later on, l/(v — z) in a power series in z/v. Therefore we 
can rewrite the Fredholm determinant if„(z _1 ) of K(z~ l ) as 

H n {z- 1 ) = # n (0)(l - iU^ 1 )), (4.7) 

wherqj 

Ruiz- 1 ) := (Q, XnHz- 1 )), Q k := ((1 - XnK(0)xn)- 1 Xn9)k (4.8) 

and Xn(k) = l[fc>n]j here the symmetry of K(Q) is being used. Accordingly 
R n (z~ r ) = (Q , Xnh{z~ x )) , as in (14.81) . decomposes as (recall that n = 2m + l) 

— 2t(z— z~ 1 ) 

Ruiz' 1 ) = Sniz- 1 ) + 6 Uz- 1 ), (4.9) 

{-z) n 

with 

S n {z~ l ) = (Q, XnHz- 1 )), UZ- 1 ) = J2 %^F- (4-10) 

fc>l ^ ' 

We set for x G Z, 

— 1 /" e -*(2-2 _1 ) 

^-^-^^T^^TT^- 1 ), 
C l( *) := ± / dz U^)Uz) (4.11) 



27Ti/ ro " (-*)*+! 



2vriiro 
C{x):=2C 1 (x) + C 2 (x). 



C2{x) := 1 M ^^ TZ^T 



s For a = (a fe ) feGZ and 6 = (6fc)fc 6 z, the inner-product (a, 6) := Y,kez akbk - 
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Remark that C\(x) = C\{—x) and C 2 (x) = C 2 (— x). Also introduce functions 
Ei(z, w), which also depend on n = 2m + 1, 

pt(Z—Z ) / -y \ra 

M*M ■= ^^Fry (-) (1 - S n (z-i))(l - S n (w)) 
E2(z,w) := ~ !lll_ 1) (-^) m (-^) m+1 (l " S^z-^Uw) 
E,{zM := -^—^(-zr^i-wy-Uz- 1 )^ - S n (w)) 

E &M ■= - eKw - w -i) {-) T n {z)T n (w-). 
With these notations, the following statement hold. 
Proposition 4.1. The kernel K. m (x,y) in i\3.13\) has the following expression 

+ 1 ' ) dz I dw St^iM ft^l + J=£_) , (4 ,3) 

as we// as the Airy kernel-like expression 
(^r- y ^^^m(^y) = C(x-y) 



A(x - c)A{y - c) + A(-x - c)A(-y - c) \ 

+ ^2 I ~ A ( x ~ c ) B (y ~ c )~ A (~ x - c) B {-y - c ) 

c>o y _b(x - c)A(y - c) - B(-x - c)A(-y - c)j 
Y^ (B(x - c)B(y - c) + B(-x - c)B(-y - c 



(4.14) 



c<0 



Proof. Let us first prove (14.131) . Consider the kernel K m (x,y) as in (13. 13ft ; 
one uses if n (z _1 ) = H n (0)(l —R n (z~ r )), as in (14. 7p . and one renames the in- 
tegration variables (w, z) — > (z, w) in the second double integral, enabling us 
to combine the two double integrals. Then, taking into account the prefactor, 



K zm ) J Fo J r e v 



w/ \ \—zr \—w 



\x+l 



x (1 - R^z-^jl - R n (w)) ^ 
z — w 
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That the single integral above equals C 2 , defined in (14. Ill) , follows from 
the fact that the —1 term can be deleted, since ^-. § v dz z y ~ x ~ l = S Xjy and 
6 Xt yl XJ t:y = 0. Multiply out (1— _R n (z _1 ))(l— -R n (u>)), use the expression (I4.10p 
of R n and the functions E^s defined in (14.121) with the result 

, . i / , / , i fi-wy- 1 (-z)y 



(2vri) 2 J ro J roz z-w\ (-z) x (-w) x +\ 

x (e x (z, w) + E 2 {z, w) + E 3 (z, w) E A {w, z)\ + C 2 (x-y). 

(4.16) 
The double integral, involving the last expression in brackets, is not in a 
usable form, in view of the saddle point method and the topology of the 
contours (see the discussion after the proof). Namely, the integrations have 
to be interchanged, at the expense of a residue term, as is given by the 
general formula (13.401) . So, using this formula, and further renaming z -H- w, 
the double integral with E± becomes 

Iff 1 f (-w) y ~ l (-z) y \ 

777-^ i dzi dw \ \ x + v ' E 4 (z, w) + 2Ci(a; - y), 

(2vn) 2 /r Jr 0z z-w\ {-z) x {-w) x+1 J 

(4.17) 

where C\{x) is defined in (14.111) . So, taking equation (I4.16P and (I4.17P into 

account, we find that formula (14.131) for the kernel K m (x,y) holds. 

Next we prove (I4.14|) . The first observation is that the kernel (I4.13P 

depends on x and y through the expression in brackets only; the latter itself 

is invariant for the interchange (x,y) h-> (— y, —x). So it suffices to consider 

the double integral associated with the first term (—w) y ~ 1 (—z)~ x only; the 

other one is automatic. Since the integration paths can be taken to satisfy 

\z\ < \w\, in the double integral of (I4.13p . one may use the series 



z — w 

c>0 



Y^ ( J , valid for \w\ > \z\, (4.18) 

—w) *-^ \ — u> ) 



and one notices that for each of the Ei, the double integral decouples into 
the product of two integrals over r : 



, dz (b dw- 
(2tti) 2 J Vq J v z-w (-z) 



r _j, e^> _ r z^(-")'-":- 2 (1 _ s , iW) 



c>0 



2vri (-z) x - m ~ cK " 7 ro 2vri e'^"^" 1 ) 

^A(x-c)A(y-c 



c>0 



(4.19) 
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To see that the second integral equals A(y — c), one performs the change of 
variable z i— ¥ 1/z. Since the only poles are at z = and z' 1 = 0, this is 
allowed; so, we do not pick up further poles. The same decoupling occurs for 
the other E^s, which yields: 

, / , E 2 (z,w)(-w) y - 1 v^ ., 



c>0 



(2vri) 2 / ro J Toz ' z-w 



(4.20) 



c>0 



and 



1 r iA *„.*<*. "OH")*- 1 



(27ri) 2 J To J Toz z-w {-z) x 
= -^B(-x + c+l)B(-y + c+l) = -^B(-x - c)B(-y-c). (4.21) 

c>0 c<0 

Then adding the same expressions with the interchange (x, y) !->■ (— y, —x) 
yields formula (j4.14p , completing the proof of Proposition 14.11 □ 

In anticipation of Section [7] on the integral representation of the limit- 
ing kernel, which will be obtained by saddle point analysis, some comments 
must be made here; they will also explain the interchange of integrals, which 
occurred in (|4.17|) . Given the future rescaling m ~ 2t with x = ^it 1 ^ 3 , 
y = ^2^ 1//3 for t — > oo, the steep descent method applied to A(x) and B(x) at 
z = —1, in particular to the part of the integrand e ±4 ^~ 2 \—z) ±m = e ±tF ^ 
respectively, uses the Taylor expansions 

F(z) :=z-z- l +2 hg(-z) = \{z + l) 3 + G(z + l) 4 , 

6 (4.22) 

log(-^) = -(z+l)--(z + l) 2 + Q(z + l) 3 . 

The steep descent path for A(x) will therefore look like v with an angle 
of approximately ±7r/3, whereas for B{x) it will look like \ with an angle 
of approximateljo ±2-7r/3 with the positive real axis. The contours of the 
four double integrals of equation (I4.13p . associated with each one of the Ei's, 
from the point of view of steepest descent analysis about z,w = — 1, are 
topologically two circles, a ^-circle inside a w-circle, which are deformed so 
that locally near z = w = — 1 they look like the set of pictures in Figure H] 
(see Section [7]), with the two circles intersecting the real axis at the common 
point z, w — — 1 and to the right of —1. 



3 The angles can be within the range 7r/3 ± 7r/6 and 2ir/3 ± 7r/6. 
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4.2 Bessel reformulation 

The purpose of this section is to express the functions A(x), B(x), C±(x) 
and C2(x), as in (14. lip in terms of Bessel functions, the expressions Q k and 
the Bessel kernel K(0), as in (14. 8 p and (14.31) . Throughout we will be using 
the integral representation of the Bessel function of order n G Z, together 
with its symmetries, 

1 r t( z — z ) 

J n (2t) = — i dz 6 — = (-l) n J_ n (2*) = (-l)»J ft (-2t). (4.23) 

27T1 J Vo Z n+1 

J n (2t) is different from the modified Bessel function I n (2t), defined in (I3.4I) . 
To do so, we shall need the following Bessel function expressions for the basic 
building blocks. 

Lemma 4.2. The kernel K(0) defined in \4-'Jty , the expressions hk and 
ge given in fl^.5| ) and \4-6\) an d the functions T n (z _1 ) and S n (z~ l ), given 
in U^TW , can be expressed in terms of Bessel functions as follows: 



K(0) k/ = Y, Ji +a+ i(4t)J k+a+1 (4t) =: B 2t (k + 1,£+ 1), g t = J m (4t) 

a>0 

-2t(z-z-^) 

hkiz- 1 ) = - J2(-zyj k+a+1 (4t) + = hiz- 1 ) + 



a>0 



e -2t(z-z- 1 ) _ ^ e -2t( 2 - 2 - 1 ) 

-,,,-r, ( — Z) k+1 



k>n ^ ' a>0 

k>n 



(4.24) 



where B t (i,j) is the Bessel kernel in l^Tj - Also, 

Q k = J2 p k,eJi+i(4t), with P k/ = ((1 - x n K^)Xn)' l )k,i. (4.25) 



On 



Proof. For K(0)k,e one uses l/(v — u) = v^ 1 J2 a>0 (u/v) a for \u\ < \v\ and 
then ( I4.23p . The same geometric series is used for h^z -1 ) in (14. 5p but with u 
replaced by z, from which the formula ( I4.24p for h k {z~ x ) and the formula for 
S n by ( ESP follow. Finally, one has g t = (-l)^ 1 J_ 1 _ £ (2t) = J e+1 (2t). D 

The more intricate term is C 2 from (14.111) . 

Lemma 4.3. The expression C 2 {x), as in ^.11\ ), equals 

C 2 (x) = l [x ^ 0] C* 2 (x), (4.26) 
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where 



C* 2 (x) = (-l) x ^- i dz-^- {Rniz- 1 ) + R n {z) - R n {z- l )R n {z)) 

^ x (4-27) 

= / „ Qk ( l[a;/0] Jk-\x\+l{4t) ~ Qk+\x\ + / , QeK{fykJ-\x\ 



k>n t>n 



Proof. One first notices that the integrand in (14 .2 7p is invariant un- 
der the mapping z i— > z^ 1 . Then, using the formula (I4.24p for 
R n (z~ l ) = YlikynQkhki.z -1 ) , one breaks up the calculation as follows. 

(a) Terms from R n (z~ v ) + R n (z). We have 

Rniz- 1 ) + R n (z) = J2 Qkihiz- 1 ) + h k (z)), (4.28) 

k>n 

and thus, by integration, one checks first for x > 0, then for x < 
and for x = 0, that, using the symmetry properties of the Bessel functions 
(see (14T23D ). 

l_j dz h k (z-^+h k (z) = { _ lY{lx>Q j k+i _ x{At) + lx<Q j k+l+x (4t)) 

= ( — l) X t[ X ^Q)Jk+l-\x\(^t)- 

(4.29) 
Substituting into the left hand side of (14 . 2 7[) gives the first term on the right 
hand side of (OTj) . 

(b) Terms from R n (z~ l )R n (z). We have 

Rniz-^R^z) = J2 QkQth k {z~ x )h t {z). (4.30) 

k,£>n 

From (147231 and (14721) we get 
1 / , h k (z- l )hi(z) 



r dz h k (z-\ht(z) = j2(-iy- b s a - b , x j k+a+ imj b+i+1 (tt) 

J V ° a,b>0 

~ ^(-l) 6+fc+1 ^ +6+1 (4t) J k+b+1+x (-At) 

b>0 

- ^(-l) a+f+1 J fc+a+1 (4t) J x - e -a-im + (-l) x S e . k , x 

a>0 

= {-l) x [b~z- k , x - ^ J k+a+ i(At)J i+a+ i-\ x \{At) J 

^ a>0 ' 

= (-l)*(fc_ M -#(0) M _ w ), 

(4.31) 
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using in the last equality the expression (14.24J) for the kernel K(0). In the 
second equality we used the symmetries (I4.23P of the Bessel functions. Sub- 
stituted into (14.301) . this gives the last two terms in (I4.27p . □ 

Proposition 4.4. The expressions A(x), B(x), C{x), defined in l[4.11\ ) for 
i6Z, can be expressed in terms of Bessel functions Jk, Qk and the kernel 
K(0), as follows: 

A(x) = J m+ i_ x (2t) + / J / ^Qfc^fc+l+q(4t) J m+ i +a - x (2t), 

k>n a>0 (4 32) 

B ( X ) = ^2QkJk-m+x(2t), 
k>n 

and 

C{x) = Y,Qk (Jk~x+i(4=t) + J k+x+1 (4t))+ Y, QkQi {K(0) k+x ,i + K{Q) k . Xtt ) . 

k>n k,£>n 

(4.33) 

Proof. The formulas for A and B follow directly from (14. lip and the expres- 
sions for T n and S n in (I4.24p . together with the symmetries f!4.23[) of the 
Bessel functions. Then 

r, is - 1 I a Tjz^Tjz) 

Ci(x) = — (p dz 



2vri J To (-ZY+ 1 

(—l) x f (— z ) k - £ 

= 2_^ QkQi r, • f dz — — = 2_^ QkQlO~k-£,x = / J QkQk+\w\- 

k,£>n ^ r ° " k,i>n k>n 

(4.34) 

From Lemma 14.31 it follows that 

C 2 (x) = 1[^ ] Yl ® k ( J k-\*\+i(*t) ~ Qk+\x\ + J2 Qi K {0)k,e-\ x \ ) • (4.35) 



k>n x i>n 



Next we show that lu^o] can actually be omitted. To do so, it suffices to 
show that the sum on the right hand side of (14.351) vanishes when x — 0. 

Indeed, setting P = (1 — XnK{fyXn)~ X , as in ( I4.25p . remember that 
gi = J^ + i(4t) and that Qk = (PXn9)k- Then, denoting (-, •) the canonical 
scalar product on £ 2 (Z) we get, for x = 0, that the r.h.s. of (14.351) equals 

(PXn9, Xn9) ~ (PXn9, XnPXnQ) + (PXn9, XnK (0)XnPXn9) 
= (PXn9, Xn9) ~ (PXn9, Xn(l ~ XnK(0)Xn)PXn9) (4.36) 

= (PXn9, Xng) ~ (PXn9, Xn9) = 0. 
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Plugging there results into C(x) = 2C\{x) + C 2 (x) we obtain 

C(x) = J2Qk[Jk-\x\+i(4t) + Qk+w + J2Qt K (°)w-\x\)- ( 4 -37) 



k>n l~>n 



It follows from the relation P = 1 + XnK (0)x„P (see the definition of Q and 
P in (I4.25P ) that acting on Xn9 and taking the fcth entry, 

Qk = l[k>n] [Jk+i(4t) + J2 K (°)wQt) • (4-38) 

^ £>n ' 

Using this relation for Qk+\ x \ in (14.371) we obtain 

C(x) = J^Qk (4-1*1+1(4*) + J fc+N+1 (4£)) 

k ~ n (4.39) 

k,i>n 

Finally, since K(0) is symmetric, we replace K(0)k t e-\ x \ = K(0)£_\ x \ t k and 
change the labeling k <-¥ £ to get (14.331) . Notice that we can then replace \x\ 
by x in this expression. □ 

5 Extended kernel for finite time 

Formula ( 14. 14j) (in Proposition 14. lj) with A(x),B(x),C(x) given by Proposi- 
tion 14.41 gives the kernel K m governing the fluctuations of the walkers near 
the point of meeting of the two groups of non-intersecting random walkers 
at time r = 0. In this section we prove Theorem 12.11 and we extend Propo- 
sition |4J] to the multitime setting (Theorem 15. 4p . 

Consider the n = 2m + 1 walks whose positions were denoted by Xk{r) in 
Section [3J Consider p different time slices T\ < t 2 < . . . < r v in the interval 
(—£,£). Then, the probability measure at these times of the positions of the 
random walks is given by 

(p n 
n riKfo) = ti}\ riKw = xk{-t) = m +i-k} 
j=i fe=i fe=i 

= const x det [p t+n (m + l-i, y))] 1 <. J < n 

•P-i 
x 



II det [Pn+t-nivl Vj +1 )] !< 4j < n ) det [p t -rM , m + 1- j)] x <. 
e=i ' ' 



3<n 

(5.1) 
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It is well known that a measure of this form has determinantal correlations 
in space-time [I21I221I2ZIEZII1S] , as stated in the following proposition. 

Theorem 5.1. Any probability measure on {x\ , 1 < i < n, 1 < £ < p] of 
the /brno 

1 P ~ 1 

— det Utro^Ti,^)) JT det \M T ti x i t) i T t+i> x ? +1) )) 

Z V J / Ki,j<n ±± V J / Ki,j<n 

xdet U^x^-.Tp+ubM , (5.2) 

\ / l<i,j<n 

has, assuming Z ^ 0, the following determinantal k-point correlation func- 
tions forti,...,t k e{n,..., r p }: 

p (fc) (ti,xi,...,t fc ,x fc ) = det(ir(t i ,x i ;t i ,x j )) 1 <. j .< fc . (5.3) 

The space-time kernel K (often called extended kernel) is given by 

K(t 1 ,xi;t 2 ,x 2 ) = -(p(t 1 ,x 1 ;t 2 ,x 2 )l(t2 > h) 

+ ^24>(ti,xi;T p+ i,b i )[B~\ j (J)(To,a j ;t 2} x 2 ) 

with (* means integration with regard to the consecutive dots) 

AU *>«■ „\- J 0(^r,a;;r r+ i,-)*---*</>(r s -i,-;r s ,y), ifr r <r s , . . 
<p[T r ,x,T„y) - I ^ %fr r >r s , [b - b) 

and with the n x n matrix B having entries Bij = 0(to, c^; Tp+i, 6j). 

Our measure (15 .ip has the form required by Theorem 15.11 The normaliza- 
tion constant Z is nothing else but the partition function and it is non-zero 
since the set of n paths satisfying the non-intersection constraint is non- 
empty. We already determined the one-time kernel for r = 0. To get the 
extended kernel one has to let the one-time kernel "evolve" by means of the 
operator of the random walk. This formulation was already present in the 
work of Prahofer and Spohn on the Airy2 process |41j . 

Lemma 5.2. The extended kernel K^f of the time- dependent point process 
fj(r, x) is given in terms of the kernel K m of the same point process at r = 
by the formula: 



K e *\t 1 ,x 1 ;t 2 ,x 2 ) = -l [t2<tl] (e^-^ n ) (x u x 2 ) + [e- hn K m e t2n j (x 1 ,x 2 ), 

(5.6) 



10 The functions (J>(t£, x; Te+i, y) themselves may in fact vary with I above. 
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where the infinitesimal generator % of the single random walk, the discrete 
Laplacian, acts on functions f as 

Hf(x) = f{x + 1) + fix - 1) - 2/(x), x e Z. (5.7) 

Comparing the first term of (j5.4p and (15. 6p . one sees a different ordering 
in the times. This is consequence of the dual transformation. 

Proof of Lemma I5.M The operator % in (15. 7p is the generator of the con- 
tinuous time process defined by the transition probability p t ix,y), in (12.11) . 
Indeed, one checks that this transition probability is given by 

p t (x,y) = e-%_,|(2£) = — <b dz = e m l(x,y) = (e m )(x,?/), 



27ri J Tq " z x -y +1 

(5.8) 
because 

dt PtK ,y) 2iri J To z*~y+ lK ' (5.9) 

= Ptix- l,y)+p t ix+l,y) -2p t ix,y) = iHp t )(x,y) 

with initial conditions p ix,y) = l(x,y). Here, 1 denotes the identity oper- 
ator on Z, i.e., l(x, y) — 1 if x = y and t{x,y) = if x ^ y. The one-point 
kernel in Section [31 formula (I3.38p . was written as a sum involving ipkix) and 
tpk(y)- Under the time flow, they will become different functions; therefore, 
we set \l/fc(0,x) = $&((), a;) = ipkix), and thus, with this new notation, the 
kernel reads 

n n 

K m (xi, x 2 ) = ^ ^fc(xi)V*(»a) = 5Z * fc (°> x i) $ fc(°' x 2)- (5.10) 

fc=i fe=i 

The two set of functions {$fc(0,x), fe = 1, . . . ,n} and {\l/fc(0,x), fe = 1, . . . ,n} 
satisfy 

span{$fc(0,x), k = 1, . . . ,n} = span{pt(m + 1 — k,x),k = 1, . . . ,n}, 

span{\I/ fc (0, x), k — 1, . . . , n} = span{p t (x, m + 1 — k), k = 1, . . . , n}, 

with ($ fc (0, x),^ p (0, x)) = 4, P , 

(5.11) 
so that the matrix B defined in (15. 4p becomes the identity matrix. 

Let us consider the functions of Theorem 15.11 First of all, the function 
4>iti,Xi,t 2 ,X2) appearing in ( 15 .4p becomes 

(j)it 1 ,x 1 - 1 t2,x 2 )t[t 2 >t 1 } = t[t 2> t 1 ]Pt 2 -t 1 ixi,x 2 ) = t[t 2> t 1 ]ie (t2 ~ tl)n )ix u X2), 

(5.12) 

29 



where ti, t 2 G {ti, . . . , r p }. Next, with r = — t, r p+1 = t we have 

0(ti, z; t, 6*) = vt-tx (x, h) = (e- tin )(x, •) * 0(0, • ; t, b k ) (5.13) 

and 

(f)(-t,a k ;t 2 ,x) =p t2+t (a k ,x) = <f)(-t, a k ;0, •) * (e* 2W )(-,x). (5.14) 

With the choice of basis used for the kernel at r = 0, we have that 0(0, ■ ; t, b k ) 
is replaced by ^(O, •) and </>(— t, a^, 0, •) by $fc(0, •) (so that B = 1). Thus 
in Theorem 15.11 we have replaced 

0(*i,s;*A) -► (e- tlW )(x,-)*^(0,-) = (e- tlW * fc (0, -))(x) =: * fc (<i,ar) 

(5.15) 
and 

<K-t,a k ;t 2 ,x) ^$ fc (0,-)*(e t2W )(-,x) = ($,(0, -)e t2W )(x) 

= (e i2wT $ fc (0,-))(x)=:<l> fc (t 2 ,x). 

(5.16) 
Therefore the extended kernel has the following expression in terms of the 
kernel K m in (137381) 

n 

K^(ti,a;i;t 2 ,x 2 ) = -l[ti<ta]Pt3-ti(«ij«2) +^**(*i,a;i)$fc(*2,E2) 

fc=i 
= -l [tl<fa] (e( t2 - tl ^)(x 1 ,x 2 ) + (e- hH K m e^ H ) ( Xl ,x 2 ). 

(5.17) 
Notice that, using the semi-group property of e tn , we have the consistency 
relations (for i — 1, . . . ,p) 

9 fax) = (e^-^* fc (r p ,-)) (s), 

$ fe (r„x)=($ fe (r 1 ,-)e (T ^ l)w )(x). 

The kernel IK^f for the dual random walk is then given by taking the 
complement. Using (I5.17P and remembering that K m = 1 — K m from for- 
mula (13.101) we get 

K^(ti, xi, t 2 , x 2 ) = l [tl=t2] l(xi, x 2 ) - KJf (ti, xi; £ 2 , x 2 ) 

=l [tl =t a] l(xi,x 2 ) + l [tl<t2] (e(* 2 - tl ^)(xi ) x 2 ) - (e" tlW K m e i2W ) (xi,x 2 ) 

=l [il=t2] ( e ( t2 - il ^)(x 1 ,x 2 ) + l [tl<t2] (e(* 2 -* 1 ^)(a;i,x 2 ) 

- (e^-^ H )( Xl ,x 2 ) + (e^ n (l - K m )e <2W ) (x h x 2 ) 

= - l [t2<tl] (e^-^ n )( Xl ,x 2 ) + (e-^ n K m e t2n )(x u x 2 ), 

(5.19) 
yielding (15. 6p . ending the proof of Lemma [5.21 □ 
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With the help of Lemma 15.21 we can easily prove Theorem 12.11 starting 
from Theorem 13. 11 



Proof of Theorem \2.1[ One of the key ingredients is that f(x) := u x is an 
eigenfunction of % with eigenvalue u + u^ 1 — 2. Indeed, 

(Hf)(x) = u x+1 + u x ~ x - 2u x = {u + u- 1 - 2)u x = (u + u- 1 - 2) f{x). (5.20) 

Moreover, 7i is symmetric. Therefore, 

(e m f)(x) = e^ u+u ^-^f(x), (fe m )(x) = (e mT f)(x) = e t{u+u '" -^ f [x) . 

(5-21) 
Then, (12. 5p follows straightforwardly from (I3.13P by applying e tlH to the 
left, e t2H to the right of K m (together with ( 15.81) for the first term of (15. 6p ). D 

For the further analysis, we extend the reformulation of the kernel for 
t = 0, as in Proposition l4.lt to the extended case. For that purpose, we first 
define the basic functions replacing A, B, and C of the one-time case (see 
Proposition |4.4p . To do so, define a new function Jx(2t) dependent on a 
parameter r 



J:Vt) := i ^——e^-- 2) = e~* ^- J X (2VW^ 



Fo 27ri2; z x \t — r 

(5.22) 
Also define a r-dependent extension of the kernel K(0)k/, as in ( I4.24p . namely 

K^(0) k/ := J2 Ja+k+Mt) Ja+e+i^t). (5.23) 

a>0 

Then define new functions A(t,x),B(t,x),C(t,x) with reR and x G Z, 
which extend the functions A(x),B(x),C(x), first defined in (14.111) and re- 
expressed in (14.321) . by 

A(r,x) := ji T |i-( 2i ) + EE^^+^^^ii+o- ( 2t )> 

fe>n a>0 

B(r,a;):=X>j£L*(20, 

fc>n 

ccr,*) := $^g fe (4 T + } 1+ ,(4t) + j&-( 4 *) 

k>n 

+ $^ Q^ (^ (r) (0) fc+ .,£ + ^ (r) (0) fc -^) • 

k,£>n 

Remember H n (Q) = det(l - K(0))^ n , n +i,...)- 
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k>n . . 

(5.24) 



Lemma 5.3. Given the notation \4-lty f or ^e Ei 's, the extended kernel K. c r 
is given by 

'-l)^tf n+1 (0)~ ext 



-l)^e 4t i H n (0) 



K^(t 1 ,x l ;t 2 ,x 2 ) 



Mt2<h]Pti-t 2 (xi, x 2 ) " + * +C(tt -t 2 ,x 1 -x 2 ) 



+ i^i dz { dw 7=zi2 %(*>") 



(( — w ) x 2-l e -t 1 (z+z- 1 +2) ( — zY 2 e - t i( w + w ~ 1 + 2 ) 

(— Z ) x i e -t2(w+w~ 1 +2) + (Z^/W+l e -t 2 (z+z- 1 +2) 

Theorem 5.4. The extended kernel K^f is also expressed as 

K e ^(t 1 ,x 1 ;t 2 ,x 2 ) 



(5.25) 



-l)^t 2 Hn+i{0) ^ 
-l)^e^ H n (0) 

= -Mt2<ti]Ph-t 2 (xi,x 2 ) " +1 +C(ti -t 2 ,X! -x 2 ) 



+ 2j f^(*ij #i - c)A(— 1 2 , x 2 - c) + A(ti, -xi - c)A(-t 2 , —x 2 - c) 

c>0 

—A(t l7 x\ - c)B(-t 2 , x 2 - c) - A(ti, —xi - c)B(-t 2 , ~x 2 - c) 
—B(ti, Xi - c)A(-t 2 , x 2 — c) — B(t%, —xi - c)A(-t 2 , —x 2 - c) 

— y^ \B{tx,x-i. - c)B(—t 2 ,x 2 - c) +B(t 1 ,—x 1 - c)B(-t 2 ,-x 2 -c)J. 

c<0 

(5.26) 
Proof of Lemma \5.3\ and Theorem\5.J\ First of all, let us focus on the term 



(eS t2 4l ^ w ) (^1,^2) i n (15. 6p . Remember that t 2 — ti< 0, so we can rewrite 

pi^7 ( e(t2 " l} ^ (—^) = pi^k^l— l(2(fa " *0) 

= ^7^|x 1 -x 2 |(2(ti-t 2 )) =p tl _t 2 (x!,x 2 ), 

(5.27) 
where we used the property I n (—2t) = (—l) n I n (2t) of the modified Bessel 
function; see (13.4J1 . 

Next we derive the double integrals in (I5.25p . The corresponding expres- 
sion of the kernel K m in (14. 13j) is a linear combination (not forgetting the 
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conjugation factor of the l.h.s. of (14.131) ) of 

(5.28) 



w X2 1 z X2 



Z x ^ 7y^l"rl 

Applying e~ tin to the left and e t2H to the right, (15.281) transforms into 

__ f . (5.29) 

z xi g-t 2 (w+-iu 1 — 2) yjxi+l e ~t2{z+z- 1 —2) ' \ ' J 

The multiplication by the prefactor r^ Xl e e i tl leads then to the expression in 
fl5T25|) . 

Next derive the terms with the sums in (15.261) and the expression for C. 
We act with the semigroup on the summation part of the kernel (I4.14p . which 
is expressed in terms of A(x), B(x), C(x), namely 

%±lMK m (x 1; x 2 ) = ^[(-1)-^ - c)][(-irA(x 2 - c)] + . . . 
H »(°> ^ (5.30) 

+(-ir- x2 c(x 1 -x 2 ), 



with A(x) , B (x) , C (x) given in Proposition 14.41 So, except for the term 
C(x\ — x 2 ), the expression above is a sum of decoupled terms. Therefore 
acting on the (— l) x A(±x — c)'s and (— l) x B(±x — c)'s with e~ tin to the left 
amounts (by linearity) to acting on the (— l) x JN± x (2t) (for some N depending 
on the terms) and finally to acting on l/(— z) ±x inside the integration. More 
precisely, by (15.211) with f(x) := l/(— z) ±x , we have 

(e- hn f)(x) = e tl{z+z ~ 1+2) f(x), and (fe t2H )(x) = e -^*+*~ x +*) f( x ) } 

(5.31) 
from which, by linearity, 



1 tlz—z^ 1 ) 
aZ 6 p t 1 {z+z- 1 +2) _ (_-\\x f> Ui j{tl) 

r 27rizz N (-z 



Y—^ ei) = (-l)^ 4tl J^(2t), (5.32) 



and 

Y,(-l) y JN± y (2t)(e^)(y,x) = (-l) x e-^J ( -H\2t). (5.33) 



This extends to the functions (— l) x A(±x — c), (— l) x B(±x — c) because they 
are linear in the (—l) x JN±x(2t) (see (I4.32p ). Explicitly, applying e~ tlH (to the 
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left) to (— l) x A(±x — c) amounts to replacing A(±x— c) with e Atl A(ti, ±x — c). 
Similarly, applying e t2H (to the right) to (— l) x A(±x— c) amounts to replacing 
A(±x — c) with e~ 4t2 A(—t2,±x — c). The same holds for B instead of A. 
Thus we have obtained the terms in kernel (15.261) including A's and B's. 

Exactly the same procedure applies for the term (— l) Xl ~ X2 C(xi — X2), 
because it is again a linear combination of (—l) Xl ~ X2 J^± xlTX2 (4:t) . Therefore 
acting with e~ tlH and e t2K as before on (— l) xl ~ X2 C{xi — x 2 ) leads to the 
replacement of C[x\ — x 2 ) by e 4 ^ 1_ ' 2 - ) C(ti — t 2 , x\ — x 2 ). This ends the proof 
of the formulas ( I5.25P and (I5.26P for the extended kernel, thus establishing 
Lemma 15.31 and Theorem 15.41 □ 



6 Asymptotics 

In this section we prove the first half of Theorem I2.2[ namely formula (I2.20p . 
From the discussion in Section [2] after Theorem 12.11 concerning the inter- 
action between the top and bottom sets of random walks, we rescale space, 
time and the gap n = 2m + 1 between the two groups of walkers, as follows: 



in 



2t + at 1 / 3 , x i = i i t 1 '\ t^Sit 2 / 3 , 2 = 1,2, (6.1) 



where a 6 R is a fixed parameter modulating the "strength of interaction" 
between the upper and lower sets of walks. To prove formula (12.201) of The- 
orem I2.2[ we first analyze the asymptotics of the building blocks and deter- 
mine some bounds which will be used later to show that we can exchange 
(by dominated convergence) the large time limit with the integrals (sums). 

Recall from ([532]) . (12~T31) and (E22D the functions J x T) {2t) and Q, and 
the kernel K {r \Q) k/ : 

j£\2t) = e"* i^j' 2 U2VW^). 
K^(0) k/ = J2 4;Vi( 4 *)<Wi(4t). ( 6 - 2 ) 

a>0 

Q{k) = [(1 - xaKMXay'XaM]^), with a := 2 2 ' 3 a, 
and where Xa{x) = l[x>a]- Remember from (I2.14p the definition of 

Ai (s) (0 :=e ?s+ F 3 Ai(£ + s 2 ), (6.3) 

and define the Airy-like kernel 

/■oo 

tf«(«,A):= / d 7 Ai (s2_2/3) ^ + 7)Ai(A + 7 ). (6.4) 

Jo 

34 



Also define the following step functions of k,A G K, for which -by 
anticipation- we indicate the limits for t — > oo: 



Jt ( K ) '--t J [2t+KtV3+l]( 2 *) 



Ai (s) ( K 



K) 



K-t ( K A) '■- (2t) 1/3 K (s * ) (0)[4t+ K (2t) 1 /3] > [4t+A(2t) 1 /3] 



'W 



Q t (/c) := (2i) 1/3 Q [4t+K(2<)1 



/31 



1 — \ n-it JCj. X n - At 
k (2t) 1 / 3 ' (2*) 1 / 3 



X^« c/2t 

(2t) 1 /3 



tftf(M) 



Q(k). 



(6.5) 



Lemma 6.1. VFe have the following bounds and limits for J7 t and /Q 
defined in A6.5\) . There exists at Q >0 such that uniformly for t >t Q it holds 
that 

\J t {s) (K)\ < Cl min{l,e- eK }, \K,\ s \k, A)| < c 2 e~ e{K+x) (6.6) 

for any fixed 9 > and some constants ci,c 2 > (independent oft). More- 
over 

lim J} s \k) = M {s \k), lim ^(k, A) = K l J$(n, A) (6.7) 



t— >oo 



t— >oo 



uniformly for k, A and s in a bounded set. 
Proof. We have 

^ ) (0=t 1 / 3 J [ ^ 1/3] (2t) 



71 ,-1/3 \ i[2i+?* 1/3 ] 



1-st- 1 / 3 
The prefactor can be estimated for t — > oo, as follows 

-l/3\'+^ tl/3 



-2st 2 / 3 



1 + St~ 

1 - St-V3 



e ^+3* J (1 + ©(r 1 / 3 )) 



(6.8) 



(6.9) 



where the (9(t _1//3 ) is uniform for s in a bounded set and independent of 
£. Therefore, for £ large enough, | (16.91) | < exp(2|£s| + |s 3 |). Concerning the 
remaining part of (16.81) . using (1A.4J) one readily obtains 



lim t 1/3 J [2mtV3] (2t V / l - * 2 £- 2 /3) = Ai(f + s 2 ). 



(6.10) 



Regarding the bound, for s in a bounded set, if t is large enough it follows 
from the bound (IA.6[) that 



i 1/3 <Wv3] (2*Vl - sH-^) 



(6.11) 
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is first of all uniformly bounded and for large £ it decays as e _/3 ^ for any 
choice of /3 > 0. The statements in the first parts of ( I6.6P and (16 .7p then 
follow if we choose (3 satisfying (3 > 9 + 2|s| for any s in the given bounded 
set. 

To compute the limit of /Q , one uses definition (16.51) and formula (I6.2p 
for if( s * )(0), but with J replaced by J in the last equality below, 

/Cj S) (M) = (2t) 1/3 K (st2/3) (0) [ 4 t+K ( 2i ) 1 /3 ] , [4i+ A(2 i )V3 ] 

= (2t) 1/3 J] J [4f+( 7 + K K2t) 1 /3]( 4t ) J [4t+( 7 +A)(2t)V3 ] (4t) 

7G(2t)- 1 /3pj (6-12) 

x; ^r 2/3) (.+7)#(A+ 7 ). 



(2t)V3 

V ' 7 G(2i)- 1 /3N 

From this, using the bound (I6.6P on J~, we obtain 

|/C^(K,A)|<c? e ^^)-i^ £ e~^<c 2 e-^ (6.13) 

^ ' 7 e(2t)- 1 /3N 

for t > to = 1 and some C2 > 0, uniformly for s in a bounded set. 

We can think of the sum in (16. 12)) as an integral of piece-wise constant 
functions. The first bound in (16. 6p allows us to use dominated convergence to 
exchange the limit and the integral. Then, lim^oo <Jt ( K ) = Ai^(«) yields 

POO 

lim K.\ s \k, A) = / rf 7 Ai (2 " 2/3s) (/€ + 7 )Ai(A + 7) = K$(k, A). (6.14) 

□ 



t— >00 JQ 



Lemma 6.2. Set a t := j^tjs an d define the operator M. t = Xat^t X&t> 
appearing in the definition Ii6.5\) of Q t . Then, uniformly for t > t , we have 
for the operator-norm n \ \\ ■ \\, 

\\Mt\\ < 1 (6.15) 

which implies that 

||(1 - M) _1 || < (1 - IIMII) -1 < C < 00 (6.16) 

for some finite constant C independent of t. 



1 where \\A\\ = sup^^ \Af\. 
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Proof. By Lemma [6. II and the fact that o t — > o as t — > oo, it follows that 

lim M t = XaK MX a =■■ M (6.17) 

t— >oo 

pointwise. Moreover, 
lim ||M-.M|| 2 < lim ||M--M|||s= lim f dndX \M t {K, A) - M(k, A)| 2 

t— >oo i— >oo t— >oo J 

= / dndX lim \M t (K, A) - M(k, A)| 2 = 
y t— >oo 

(6.18) 

where we use by Lemma 16.11 dominated convergence to exchange the limit 

and the integral together with (J6.17)) . It is known that A max = ||A^|| < 1 for 

any fixed a (see, e.g., [IH]). This, together with (j6.18p . implies that 

\\M t \\ < \\M\\ + \\M t - M\\ < 1 (6.19) 

for t large enough. D 

Lemma 6.3. Consider Q t as defined in Ii6.5\) . There exists a t > such 
that uniformly for t >t it holds 

\Qt(K)\<c 3 e- dK (6.20) 

for any 9 > and some constant C3 > (independent oft). Moreover, 

lim Q t («) = Q(k) (6.21) 

t— >oo 

uniformly for k in a bounded set. 

Proof. For the sake of this proof, set J t := Jt an d ^Q := ^Q ■ First 
of all we prove that Qt{n) is uniformly bounded for t > to- Recall that 
<2 t (/t) = [(1 — A^i) _1 X5- t J2t]( K )- Since (1 — M. t )~ l exists, we can use the 
identity 

(1 - Mt)' 1 = 1 + x*,£tx*t(l - M) -1 , (6.22) 

which upon integrating from a to 00 against the function J 2 t gives 

/■oo 

Qt(«)=X*t^(«)+ / d\K t {K,\)\{\- M t y\~ at 3^{\). (6.23) 

Thus, 

^■00 
|Qt(«)|<|x*A(«)|+ / rfA|/C t (K,A)||[(l-M i ) _1 X^^](A)|. (6.24) 

37 



But 

|[(1 - M t )- l X , t J 2t }{\)\ < ||(1 - M t y l \\ 1^4 (6.25) 

is uniformly bounded for t > to (by Lemma 16.11 and Lemma 16.21) . Then, 
using the bound for K, t and J t (k) in (16. 6p we obtain the bound (16.201) . 
To prove (16.211) . we show that 

I Qt ~ Q\oo = sup | Q t («) - Q(k) I -> 0, (6.26) 

as £ — )• oo. We have 

IQt - Qloo = |(1 - M t )- l XiJt J 2t - (1 - A*)- 1 **^ 

< | [(1 - M t y l - (1 - A^)- 1 ] x.^tL (6-27) 

+ |(i - m)- 1 ixaJv - x*ti\L + <?r l/3 ), 

where the correction term 0(t~ 1 ^) comes from the fact that the difference 
between a t and a is not larger than (2t)~ 1//3 . Then, 

dS2ZD < ||(1 - Mt)- 1 - (1 - A^)- 1 !! \xaJ2tL 

+ ||(i-A^)- 1 || \ X9 J* - x^L + 0(r 1/3 ) J '" 



The first term goes to zero as t — > oo. Indeed, \XaJ2t\oo < C < oo by 
Lemma 16.11 and, using the identity 

(1 - Mt)' 1 - (1 - MY 1 = (1 - Alt)" 1 [Al t - Al] (1 - X)- 1 (6.29) 

together with the fact that \\M t \\ < 1, ||Af|| < 1, and ||Af-Af t || ->■ 
in the t — >■ oo limit (see Lemma 16.21 and (I6.18P ); so one has 
||(1 — Aft) -1 — (1 — Al) _1 || —?• 0. The second term goes to zero as well, since 
||(1 — Af) -1 || is bounded and, by Lemma ISTTj \xdJ2t — Xa-Ai|oo ~~ ** 0- D 

Proof of Theorem \2.2\ formula l[2.2U\) . We now define new functions 
At(s,£), Bt(s,£), Ct(s,£), which are rescaled versions of A(r, x), B(t,x), 
C(t,x) (see formula (I5.24p ) under the scaling ( 16. ip : 

A(s,0 



t 1 / 3 A(st 2 / 3 ,^ 1/3 ), 

t^B^' 3 ,^ 3 ), (6.30) 

t^ 3 C(st 2 / 3 ,^ 3 ). 



As t — > 00, these functions will converge to A(s, £),B(s, £),C(s, £) of (I2.15P 
and (I2TT6D . 
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One then recognizes in these expressions the functions (16.51) . thus yielding 
+ 72*173 E^ E QMjS\- + ^j} S \^a + a-0, 



KG/, 



n , / 



2 -l/3 



^.o = 4w3 E fit(«)(^ 2 " 2/3a) («- 2- 1/3 o +jT /3s W 2-^0 



KGln.t 



2-1/3 

+ 



^ £ Q t (^)Qt(A)(/C t (s) ( / ,-2- 1 /3 e)A ) + /c! s) ( K + 2- 1 / 3 e,A)). 

(6.31) 
For instance, the function Jfc + i +a (4t) in A(r, x) becomes, upon setting 
a = «(2t) 1 / 3 and « := (2t)- 1 / 3 (A; - 4£), 

J fc+ i+a(4t) = J [4i+ (« +Q )(2t)V3 +1] (4t) = (2t)-V 3 jg\K + a). (6.32) 

Notice that the sum over k > n in the expressions (I5.24p becomes a sum over 
« £ -fn.* with 

J ftit := (2t)- 1 / 3 ({n, n + 1, . . .} - it). (6.33) 

so that the condition k > n = 2m + 1 = At + 2er£ 1 / 3 + 1 translates into 
K = (2t)~ 1 ^ 3 (k — At) > 2 2 / 3 ct = a. Setting the summation variable c = 7t 1 ^ 3 , 
rewrite the kernel (I5.26P in Theorem I5.4[ with the scaling (16. ip 

(-l)*>e^ H n+1 (0) ~ ext 



;-i)^ e ^ fl- n (o) 



K^ t (i 1 ,x 1 ;i 2 ,x 2 ) 



/ A(si,6-7)A(-s 2 ,6-7) + A(si,-ei-7)A(-s 2 ,-6-7)\ 

+ 7T/J E -A(si,6-7)^(-s 2 ,6-7)-A(si,-6-7)^(-s 2 ,-6-7) 

7et -i/s N y - B t ( Sl , 6 - 7 ) A(-«2, 6 - 7) - Bt( Sl , -6 - 7) A(-s 2 , -6 - 7)/ 

- 7771 E { B ^ & ~ 7)fi*(-*2, 6 - 7) + BtK -6 - 7)Bt(-«2, -6 - 7)) 

7Gt" 1 /3z_ 

(6.34) 
In view of (I2.10p we have lim t ^ oo if n+1 (0)/if n (0) = 1 and in the t — > oo 
limit, (n — 4t)/(2t) 1,/3 — > a. Notice that the sums with the preceding volume 
element, 1/t 1 / 3 or l/(2t) 1 / 3 depending on the case, can be just thought of as 
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integrals with the integrand being piece- wise constant. What follows holds 
uniformly in t for t > t where to is a fixed constant. The exponential bounds 
of Lemma 16.11 and Lemma 16.31 imply that for any 9 > there exists some 
c > (the constant c depends on a, which is however fixed) 

\A t (s,-£)\ <ce~^ and lim A t (s,£) = A(s,£). (6.35) 

t— >oo 

Moreover At(s,£) tends to A(s,£) uniformly on bounded sets, by uniform 
convergence on bounded sets and dominated convergence of the integrand. 
Using the exponential bound of Lemma 16.31 and the fact that Jt is just 
bounded, we obtain similarly 

|B t (s,OI <cmin{l,e - ^} and lim B t (s,£) = B{s,$,). (6.36) 

t— >oo 

Finally, the exponential bounds of Lemma 16.11 and Lemma 16.31 imply that 

\C t (s,0\<c and lim C t (s,0 = C(s,0, (6.37) 

t— >oo 

where the last limit holds uniformly for £ and s in bounded sets. 

Using the bounds in (16.351) . f)6.36p . and (16.371) . one concludes that the 
integrands (summands) in (I6.34p are uniformly bounded by functions which 
are integrable (summable). This is uniform for £,77 and s in a bounded 
set. Then, by dominated convergence, we can take the limit inside, thus 
yielding (12.201) . Finally, the Gaussian term in ( I2.20p comes from the known 
asymptotic (for s > 0), 

lim t l ^e- 2st2/3 L t , /3 {2st 2 ^) = ^exp K7(4*)) , (6.38) 

which can be derived from a saddle point argument. □ 

7 Integral representation of the Tacnode ker- 
nel 

To derive the double integral representation ( I2.2ip of Theorem 12.21 there are 
two ways. One can use the Airy functions integral representations ( 1A.7|) 
together with 

f°° 1 

/ d\e~ x{u - v) = , whenever U(u - v) > 0. (7.1) 

Jo u-v 

This is quite straightforward but it requires several computations which are 
not reported here. 
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The second is to do a steep descent analysis starting from formula ( I5.25P in 
Lemma 12.11 Here we merely indicate a sketch of the saddle point argument 
(not a proof). The limits of the other terms have been discussed in the 
previous section. The main task here is to take the limit of this double 
integral, when t — y oo, with the scaling 

n = 2m+l, m = 2t + at 1 ' 3 

z = -l + (r 1/3 and w = -l + wr 1/3 (7.2) 



■■*■- y'. 



&t 1/3 and t i = s i t 2/3 , z = l,2. 



Also recall the definitions (12.17P of the Laplace transforms <2(C) and V((), 
as well as the function C in f)2.16p . The reader is reminded of the steepest 
descent discussion in Section fl~Tl For taking the limit of the extended kernel, 
we need the following Lemma. 



Lemma 7.1. Given the scaling \l.ty above, the following limits hold: 

lim e t{z - z ~ X \-z) m = e^"* (7.3) 



t— >oo 



and 



lim T^z- 1 ) = e- 2 ^Q(C), lim T n (w) = e 2 ™Q(-u) 

t— >oo t— >oo (7 A\ 

lim S^z- 1 ) = V(0, lim S n {w) = V(-u), 

t— KX3 i— >oo 



where V and Q are the Laplace transforms defined in ^2.11\). One also checks: 



lim ( . W) * 2 - = ^ «nd lim e^ (z+ ^ 1+2) = e^ 2 . (7.5) 

t^oo ( — z) Xl e^ 2UJ t^oo 

Proof. Letting t — Y oo; setting n = 2m + 1, m = 2t + at 1 / 3 , the critical point 
will be at z, w — —1 and thus the leading contribution will come from the 
neighborhood of the critical points, which suggests the scalings in z and w 
above. The Taylor expansion of the F- function (I4.22p gives 

e t{z-z~ x ) ( z \m _ e t{z-z- 1 )+rnlog(-z) _ e tF{z)+at 1 / 3 log(-z) 

= e *F(-i+Ci- 1/3 )+^ 1/3 iog(i-Ci- 1/3 ) = e^t-^Cfi + 0(r 1/3 )). 

(7.6) 
Setting in addition the scaling for ti and Xi, one finds by Taylor expanding 
about z = —1 and w = —1 the limits ( 17. 5p . Introducing the running variable 
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k = 4£ + k(2£) 1 / 3 , one gets 
lim TJz- 1 ) = lim V - Q u k - = lim V Q k e~ ik - n+1 ^° s ^ 

t^oo V ; t^oo^ - z )k-n+l t^ooA^^K 

k>n v ' k>n 

= lim (2tr 1/3 £ (2i) 1/3 Q 4t+K(2t )V3e-^^ 1/3 (-^ 1/3 ) 

K>CT+(2t)" 1 /3 

= /" d«:Q(«:)e( K -^ 2l/3 =e- 2 ^Q(C) 

J K>a 

(7.7) 
and similarly 

lim T n (w) = e 2auJ [ dn Q{K)e- KU)2l/3 = e 2aw Q(-co). (7.8) 



t— >oo 

J K>a 



The limit of the expression S n , as in (I4.10p . involves hk, as in (I4.24p . Using 
the formula ( I4.24p for /^(z -1 ) in terms of Bessel functions and Lemma [6.11 
one checks, introducing the running variable a = /i(2t) 1 / 3 , 



lim h k {z x ) = - lim V {-z) a J k+a+1 {4t) 

t->oo t— >oo *-—' 

a>0 

= - lim (2t)- x / 3 V e^ 1/S W-Or^jWfa + ,A 

K>5-+(2t)" 1 /3 
/•oo 

= - / dfie-^ 2l/3 Ai(K + fi). 
Jo 

Therefore, one finds 

lim Sniz' 1 ) = lim (Q^nHz^ 1 )) = lim VQt/it^ 1 ) 



(7.9) 



fc>r 

hm(2t)- 1 / 3 V (2t) 1 / 3 Q t (^)h k (z- 1 ) 

t—too *— ' 

K>CT+(2t)~ 1 /3 

-/ d«Q(«) / rf/ie- K2l/3 Ai( K + /i) =P(C). 
Jk>5- Jo 



(7.10) 



This ends the proof of Lemma 17.11 □ 

Sketch of Proof of Theorem \2.2l formula Ii2.21\) . Since the sum in brackets 
in (J5.25P is invariant under the involution X\ <r> —x<i and t\ -H- — ^2, it suffices 
to consider the double integral, with the first term only. The second half 
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i / z w 




Ei(z,w) ~ 

t(F(z)-F(w)) 



E 2 (z } w) ~ 

p t(F(z)+F(w)) 



E 3 (z,w) ~ 

e -t(F(z)+F(w)) 




t(F(z)-F(w)) 



Figure 4: Contours z6To and w G To yZ in the neighborhood of z 
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Figure 5: Vertical lines ±5 + iR and ±25 + iR of integration for £ and a;. 



comes for free by acting with the involution! Given scaling (17. 2p . Lemma I7TH 
yields 



lim 

t— >oo 



1 / 3 dz dw (—w) 



a=2-l p ~ti(z+z- 1 +2) 



dC, doj 



z-w (~z) Xl e - t2( - w+w 1+2 ) ( - uj \e^ J e S2UJ 



MC 



3 «iC 2 



(7-11) 



and from (H~T2~j) 

4 

lim y^Ei(z, 



1=1 



■fl> 



^ 3 ^ 



1-P(0)(1 -£(-«)) 






e 3 



e 3 



+*C 



-e- 2 ^(l-^(-o;))Q(C)- 



e 3 






e 3 



e 2 ^{l-V{Q)Q{-uj) 

2CT .2(-0QM- 

(7.12) 

Combining (17. lip and (I7.12p yields the following limit below, first with 
the contours as indicated in Figure HJ which then can be transformed into the 
vertical lines above in Figure |5j compatible with Figure HI Indeed, to pick 
steep descent paths about z = w = — 1 respecting the integration contours 
in <f r dz § v dw of (I7.13p . one must choose the local paths, as illustrated in 
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Figure HI these paths must be completed by closed contours encircling the 
origin deformed to provide steep descent contours. In the (, uo scale, there are 
4 rays emanating from the origin uo = ( = 0; one is then free to deform these 
rays so as to obtain two parallel imaginary lines near the origin, as depicted 
above. Therefore the following limit holds for the first double integral 

lim — (t> dz (b ; r — — ; :— r > EAz,W) 

t^oo (27ri) 2 7 ro J Tqz z-w (-z) x ^ e - t ^ w + w ' 1 + 2 ^ ^ 

1 f dC f duJ gjll^l fe* lC \ (l-V(Q)(l-V(-u)) 

** I til** _ Co/,)^ 



(27ri) 2 J s+ .^ J_ s+iM e f-™ e^ 2 \e&*J t-u 



£t-~f .^a 



(2vri)2 7 2(5+iR s J s+m e -f-au> e^ 2 \e^J ( - uo 



.si 



_ _i_ / d< / to 4^i£!l (41) %OSH . {iv) 

(2jti) 2 A +11 s J- i+1 « e N?+-e«" 3 Ve fe V C"W 

(7.13) 
In view of the scaling f lT.2|) . the involution xi O —£2 and ii -H- —£2 induces 
the involution ^ <->■ — £ 2 and Si -H- — S2, so that the limit of the other double 
integral is given by the same formula (17.131) above, but with 

6 ** ~6 and si •(->• -s 2 . (7.14) 

We are also allowed to interchange the integration variables ( -H- — a;, 
provided the contours of integration are modified accordingly; this last inter- 
change implies 



dC, j duo remains (7-15) 

8+iR J-8+m. 

/ dQ I duo and d( duo interchange. (7-16) 

J28+m Js+M J-5+iM. J-28+iR 

So, the three combined maps, 

C <H> -00, 8i *+ -8 2 , 6 ^ -6, (7-17) 

have the following effect on the four double integrals (i), . . . , (w>) in ()7.13p : 
double integral (i) with ^7 — > same double integral (i), except for e _^ 
double integral (ii) with ^7 — >■ same double integral (iii), except for e _ s ^ 
double integral (iii) with ^7 — y same double integral (ii), except for e _ e ^ 
double integral (iv) with ^r — > same double integral (iv), except for e _ s ^ . 
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Therefore the limit 



t 1/3 1,1 dw A „ , 
urn - r— <p dz (b > EAz,w) 



(C_ u ,^2-l e -t 1 (z+z- 1 +2) ( — z) X2 e - t i{ w + w ~ 1 +' 2 )\ 

(—z) x i e ~t 2 (w+w- 1 +2) + (_ w y.i + l e -t 2 (z+z- l +2) I l'- 18 ) 

is given by r.h.s. of (J7.13P with the replacement 

+ TTw (7-19) 



Finally, in order to change the sign of the last integral, one switches the sign 
u — > —u and ( — > — (, which changes 

d( / du- into + / d( / du- . (7.20) 

Renaming variables ( — > u, oj — > v gives formula (12.21(1 . 



A Some properties of Bessel and Airy func- 
tions 

Let us recall that the Bessel function representation of order ngZ 

i r ptiz-z- 1 ) 
J n (2t) = —& dz——- (A.l) 



2?ri 7 ro 



z 



ra+1 



has the symmetries 

J n (2t) = (-l) n J_ n (2t) = (-l)"J n (-2t). (A.2) 

Moreover, 

, , j, ^,-,^-, = n_^Y Jn ( 2v « . (A . 3) 

2iri J Fo z z n \b — a/\ J 

It is well-known [1] that 

hmt^j (2t ) =A i(0. (A.4) 

t— »oo L ' 
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An uniform bound obtained in [34J is 

\{2t) 1/3 J n {2t)\ <c, c = 0.785..., neZ. (A.5) 

This bound, together with uniform expansion which can be found in [TJ is 
used in Lemma A.l of [23J to get the following result. Fix any 9 > 0. Then, 
there exists a constant to > and a constant C > such that, uniformly in 

t>to, 

|i 1/3 <Wi/ 3 ](2£)| < Cmin{l,e-^}. (A.6) 

Actually, the statement of Lemma A.l of [23] is for 9 = 1/2 but inspecting 
the proof it is straightforward to see that it holds for any fixed 9 > 0. The 
Airy function has, between others, the following two integral representations. 
For any 5 > 0, it holds 

Ai(x) = — f due u3/3 ~ ux , Ai(ar) = — / dve~ v3/3+vx . (A.7) 
27n Js+wl 2vn J-6+m 

Moreover, for any 5 > 0, it holds 

AiW(x) = e sa;+2s3/3 Ai(x + s 2 ) = — [ due u3/3+u2s - ux , 

2tti y 5+iR 

AiW(x) = e sz+2s3/3 Ai(x + s 2 ) = — / dve- v3/3Ws+vx . 

2rn J_ s+jM 



(A.8) 



References 

[1] M. Abramowitz and LA. Stegun, Pocketbook of mathematical functions, 
Verlag Harri Deutsch, Thun- Frankfurt am Main, 1984. 

[2] M. Adler, M. Cafasso, and P. van Moerbeke, From the Pearcey to the 
Airy process, arXiv:1009.0683. 

[3] M. Adler, J. Delepine, and P. van Moerbeke, Dyson's nonintersecting 
brownian motions with a few outliers, Comm. Pure Appl. Math. 62 
(2010), 334-395. 

[4] M. Adler, P.L. Ferrari, and P. van Moerbeke, Airy processes with wan- 
derers and new universality classes, Ann. Probab. 38 (2008), 714-769. 

[5] M. Adler, K. Johansson, and P. van Moerbeke, In progress. 

[6] M. Adler, N. Orantin, and P. van Moerbeke, Universality for the Pearcey 
process, Physica D 239 (2010), 924-941. 

46 



[7] M. Adler and P. van Moerbeke, String- orthogonal polynomials, string 
equations and 2-Toda symmetries, Comm. Pure Appl. Math. 50 (1997), 
241-290. 

[8] M. Adler and P. van Moerbeke, Integrals over classical groups, random 
permutations, Toda and Toeplitz lattices, Comm. Pure Appl. Math. 54 
(2000), 153-205. 

[9] J. Baik, G. Ben Arous, and S. Peche, Phase transition of the largest 
eigenvalue for non-null complex sample covariance matrices, Ann. 
Probab. 33 (2006), 1643-1697. 

[10] J. Baik, P.A. Deift, and K. Johansson, On the distribution of the length 
of the longest increasing subsequence of random permutations, J. Amer. 
Math. Soc. 12 (1999), 1119-1178. 

[11] J. Baik and E.M. Rains, Limiting distributions for a polynuclear growth 
model with external sources, J. Stat. Phys. 100 (2000), 523-542. 

[12] J. Ben Hough, M. Krishnapur, Y. Peres, and B. Virag, Determinantal 
processes and independence, Probability Surveys 3 (2006), 206-229. 

[13] A. Borodin, Biorthogonal ensembles, Nucl. Phys. B 536 (1999), 704-732. 

[14] A. Borodin and M. Duits, Limits of determinantal processes near a tac- 
node, arXiv:0911.1980; To appear in Ann. Inst. H. Poincare (2010). 

[15] A. Borodin and J. Kuan, Random surface growth with a wall and 
Plancherel measures for 0(oo), Comm. Pure Appl. Math. 63 (2010), 
831-894. 

[16] A. Borodin and A. Okounkov, A Fredholm determinant formula for 
Toeplitz determinants, Integr. Equ. Oper. Theory 37 (2000), 386-396. 

[17] A. Borodin, A. Okounkov, and G. Olshanski, On asymptotics of 
Plancherel measures for symmetric groups, J. Amer. Math. Soc. 13 
(2000), 481-515. 

[18] A. Borodin and S. Peche, Airy kernel with two sets of parameters in di- 
rected percolation and random matrix theory, J. Stat. Phys. 132 (2008), 
275-290. 

[19] A. Borodin and E.M. Rains, Eynard-Mehta theorem, Schur process, and 
their Pfaffian analogs, J. Stat. Phys. 121 (2006), 291-317. 



47 



[20] S. Delvaux and A. Kuijlaars, A phase transition for non-intersecting 
Brownian motions, and the Painleve II equation, arXiv:0809.1000 
(2008). 

[21] S. Delvaux, A. Kuijlaars, and L. Zhang, Critical behavior of non- 
intersecting Brownian motions at a tacnode, arXiv: 1009.2457 (2010). 

[22] B. Eynard and M.L. Mehta, Matrices coupled in a chain. I. Eigenvalue 
correlations, J. Phys. A 31 (1998), 4449-4456. 

[23] P.L. Ferrari, Polynuclear growth on a flat substrate and edge scaling of 
GOE eigenvalues, Comm. Math. Phys. 252 (2004), 77-109. 

[24] P.L. Ferrari and H. Spohn, Step fluctations for a faceted crystal, J. Stat. 
Phys. 113 (2003), 1-46. 

[25] T. Imamura and T. Sasamoto, Fluctuations of the one- dimensional 
polynuclear growth model with external sources, Nucl. Phys. B 699 
(2004), 503-544. 

[26] K. Johansson, Non-intersecting paths, random tilings and random ma- 
trices, Probab. Theory Related Fields 123 (2002), 225-280. 

[27] K. Johansson, Discrete polynuclear growth and determinantal processes, 
Comm. Math. Phys. 242 (2003), 277-329. 

[28] K. Johansson, Non-intersecting, simple, symmetric random walks and 
the extended Hahn kernel, Ann. Inst. Fourier 55 (2005), 2129-2145. 

[29] K. Johansson, Random matrices and determinantal processes, Mathe- 
matical Statistical Physics, Session LXXXIII: Lecture Notes of the Les 
Houches Summer School 2005 (A. Bovier, F. Dunlop, A. van Enter, 
F. den Hollander, and J. Dalibard, eds.), Elsevier Science, 2006, pp. 1- 
56. 

[30] K. Johansson, Two groups of non-colliding Brownian motions, In prepa- 
ration (2011). 

[31] S. Karlin and L. McGregor, Coincidence probabilities, Pacific J. 9 (1959), 
1141-1164. 

[32] M. Katori and H. Tanemura, Symmetry of matrix-valued stochastic pro- 
cesses and noncolliding diffusion particle systems, J. Math. Phys. 45 
(2004), 3058-3085. 



18 



[33] M. Katori and H. Tanemura, Noncolliding Brownian Motion and Deter- 
minantal Processes, J. Stat. Phys. 129 (2007), 1233-1277. 

[34] L.J. Landau, Bessel functions: monotonicity and bounds, J. London 
Math. Soc. 61 (2000), 197-215. 

[35] R. Lyons, Determinantal probability measures, Publ. Math. Inst. Hautes 
Etudes Sci. 98 (2003), 167-212. 

[36] T. Nagao, Dynamical Correlations for Vicious Random Walk with a 
Wall, Nucl. Phys. B 658 (2003), 373-396. 

[37] T. Nagao and P.J. Forrester, Multilevel dynamical correlation functions 
for Dysons Brownian motion model of random matrices, Phys. Lett. A 
247 (1998), 42-46. 

[38] T. Nagao, M. Katori, and H. Tanemura, Dynamical correlations among 
vicious random walkers, Phys. Lett. A 307 (2003), 29-33. 

[39] A. Okounkov and N. Reshetikhin, Correlation function of Schur process 
with application to local geometry of a random 3-dimensional Young 
diagram, J. Amer. Math. Soc. 16 (2003), 581-603. 

[40] A. Okounkov and N. Reshetikhin, Random skew plane partitions and 
the Pearcey process, Comm. Math. Phys. 269 (2007), 571-609. 

[41] M. Prahofer and H. Spohn, Scale invariance of the PNG droplet and the 
Airy process, J. Stat. Phys. 108 (2002), 1071-1106. 

[42] T. Sasamoto and T. Imamura, Fluctuations of a one- dimensional 
polynuclear growth model in a half space, J. Stat. Phys. 115 (2004), 
749-803. 

[43] B. Simon, Orthogonal Polynomials on the unit circle. Part 1: Classical 
Theory, Coll. Publ., vol. 54, American Math. Society, 2004. 

[44] A.B. Soshnikov, Determinantal random fields, Encyclopedia of Mathe- 
matical Physics (J.-P. Francoise, G. Naber, and T. S. Tsun, eds.), Else- 
vier, Oxford, 2006, pp. 47-53. 

[45] H. Spohn, Exact solutions for KPZ-type growth processes, random ma- 
trices, and equilibrium shapes of crystals, Physica A 369 (2006), 71-99. 

[46] G. Szego, Orthogonal polynomials, 3th ed., American Mathematical So- 
ciety Providence, Rhode Island, 1967. 

49 



[47] C. Tracy and H. Widom, The Pearcey Process, Comm. Math. Phys. 263 
(2006), 381-400. 

[48] C.A. Tracy and H. Widom, Level-spacing distributions and the Airy 
kernel, Comm. Math. Phys. 159 (1994), 151-174. 

[49] C.A. Tracy and H. Widom, Correlation functions, cluster functions, 
and spacing distributions for random matrices, J. Stat. Phys. 92 (1998), 
809-835. 

[50] C.A. Tracy and H. Widom, Nonintersecting Brownian Excursions, Ann. 
Appl. Prob. 17 (2007), 953-979. 

[51] P. van Moerbeke, Random and Integrable Models in Mathematics and 
Physics, CRM-lectures, Springer, 2010, pp. 1-130. 



50 



